UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


SRI  International 


6c.  ADDRESS  (City.  Sn re,  end  ZIP  Code! 
'33  Ravenswood  Ave . 
Menlo  Park,  CA  94025 


DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited 


5  MONITORING  ORGANIZATION  REPORT  NUMBERIS) 

AFOSH.tr.  9  0  -  0  544 


ON  1 6b.  OFFICE  SYMBOL  I  7a  NAME  OF  MONITORING  ORGANIZATION 


Air  Force  Office  of  Scientific  Research 
USAF.  AFOSR 


7b.  ADDRESS  ICity,  Stare,  end  ZIP  Code I 

Bldg.  410 

Bolling  AFB ,  D.C.  20332-6448 


8a.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

Air  Force  Office  of  Scientific. 


6b  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
lit  epphceblo) 

AFOSR/NP  F49620-85-K-Q013 


10.  SOURCE  OF  FUNDING  NUMBERS 


8c  ADDRESS  (City,  Stere,  end  ZIP  Code I 

Bldg.  410 

Bolling  AFB,  D.C.  20332-6448 


11  TITLE  t Include  Security  Oeailicetionl 

PLASMA  CLOAKING:  AIR  CHEMISTRY,  BROADBAND  ABSORPTION,  AND  PLASMA  GENERATION 


12.  PERSONAL  AUTHORISI 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO 

ACCESSION  NO. 

0  llO*  f 

*30  i 

A* 

13a.  TYPE  OF  REPORT 

Final 


16.  SUPPLEMENTARY  NOTATION 

SRI  Project  3656 


17  COSATI  CODES 


13b.  TIME  COVERED 
FROM5/15/85  T 


14  DATE  OF  REPORT  (Yeer,  Month.  Oey) 

1990  February 


15  PAGE  COUNT 


18  SUBJECT  TERMS  fConrmu#  on  r*v*rg§  if  nocetury  find  identify  by  block  number! 

^Atmospheric  plasma,  broadband  absorber,  cloaking  device, 
Epstein  profile,  noble  gases,  plasma  generation*  -ftftff" 


GROUP 


20  09 


17  0 


lYjj  ABSTRACT  I Continue  on  reverse  if  necessary  end  identify  by  block  number) 

v  Plasmas  with  electron  number  densities  from  1^  cm‘,3  to  5  *  10{^  cm '^generated  in  air  or  a  nobie  gas 
constitute  tenuous  plasmas.  If  the  ambient  pressure  is  near  atmospheric,  so  that  the  gas  density  is  of  the 
order  of  10j6  cm'^  to  1019  cm'3,  then  the  plasma  is  coilisicnal,  and  the  cold  coilisional  dispersion  relation 
predicts  electromagnetic  characteristics.  A  plasma  generated  by  an  electron  beam  or  by  photoionization 
produces  a  spatially  decreasing  electron  density.  This  spatial  electron  density,  approximated  ac  an  Epstein 
profile,  results  in  low  backscatter  and  high  attenuation.  Bandwidth  and  attenuation  are  quantified  An  air- 
chemistry  deionization  solution  is  described  that  estimates  plasma  lifetime.  The  lifetime  and  electron  density 
axe  sufficient  to  estimate  the  power  required  by  a  plasma  absorber.  Plasmas  generated  in  air  require  a  high 
input  power.  A  plasma  generated  in  a  noble  gas,  seeded  with  a  readily  ionized  material,  and  confined  by  a 
membrane,  requires  much  less  power,  and  is  technically  feasible.  '  Plasma  absorber  would  exceed  the 
capabilities  of  carbon— based  absorber  by  40  dB  and  is  broadband  from  VHF  through  X— band.  This  plasma 

can  be  used  as  a  switchable  cloaking  device  to  replace  carbon— based  absorber. 


20  DISTRIBUTION  AVAILABILITY  OF  ABSTRACT 
(3  UNCLASS  IF  IED/UNLIMITEDp0  SAME  AS  RPT 


22i  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Robert  Barker 


DO  FORM  1473,  samar  83  aps 


21  ABSTRACT  SECUI-:  TY  CLASSIFICATION 

□  otic  users!  UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code!  22c  OFFICE  SYMBOL 

(202)  767-5011  AFOSR/NP 


n  m.y  b«  U«M  unt:t  »»n»Uit«a  SECURITY  CLASSIFICATION  OF  THIS  PAGE 


83  APR  *cJition  m»y  b«  u9*a  until  »xh*uiTta 
All  othtr  »ditioni  §r«  ct>bOi«t« 


UNCLASSIFIED 


I  Mi 


Final  Report  •  February  1930 


PLASMA  CLOAKING:  AIR  CHEMISTRY, 
BROADBAND  ABSORPTION,  AND  PLASMA 
GENERATION 


Robert  J.  Vidmar,  Senior  Research  Physicist 
Remote  Measurements  Laboratory 


Prepared  for: 

Air  Force  Office  of  Scientific  Research 
Directorate  of  Physical  and  Geophysical  Sciences 
AFOSR/NP,  Building  410 
Boiling  AFB,  Washington,  D.C.  20332-6448 


Contract  F49S20-85-K-0013 
SRI  Project  8656 


Approved  for  public  release;  distribution  unlimited 


Approved  by: 

TAYLOR  W.  WASHBURN,  Director 
Remote  Measurements  Laboratory 

JOSEPH  J.  EASH,  Vice  President 
System  Technology  Division 


333  Ravenswood  Avenue  •  Menlo  Park,  CA  94025-3493  •  (415)326-6200  •  FAX:  (415)  326-5512  °  Teiex  334486 


CONTENTS 


DD  FORM  1473  i 

LIST  OF  TABLES .  v 

I  INTRODUCTION .  1 

II  TECHNICAL  ISSUES  . .  3 

A.  Air  Chemistry .  3 

B.  Broadband  Absorption .  4 

C.  Plasma  Generation  . .  5 

III  CONCLUSIONS  .  7 

APPENDIX  A  ATMOSPHERIC  REACTION  RATES .  9 

APPENDIX  B  PLASMA  LIFETIME  DATA .  29 

APPENDIX  C  AIR  CHEMISTRY  CODE  ORGANIZATION  ...  35 

REFERENCES .  43 

ATTACHMENTS 


"On  the  Use  of  Atmospheric  Pressure  Plasmas  as  Electromagnetic 
Reflectors  and  Absorbers"  (Reference  14) 

"Inflatable  Target  Support,  for  RCS  Measurement"  (Reference  26) 
Air  Chemistry  Code  Software  (inside  back  cover) 


ooesslon  For 


RTIS  GRAM 
OTIC  TAB 
Unannounced 
Justif icatlon 


By - 

Distribution/ 


Ava liability  Codes^ 
Avail  end/or 
Die*  !  Special 


TABLES 


A-l  Species  in  Simulation  Code .  9 

B-l  Plasma  Lifetime  in  Helium  with  100  ppm  Air  .  29 

B-2  Plasma  Lifetime  in  Helium  with  10  ppm  Air .  30 

B  -3  Plasma  Lifetime  in  Helium  with  1  ppm  Air .  31 

B-4  Plasma  Lifetime  in  Air  .  32 

C— 1  Program  Organization .  36 

C-2  Program  Elements .  38 

C-3  Auxiliary  Software  .  40 


v 


I  INTRODUCTION 


This  report  summarizes  the  results  of  a  three-year  investigation 
conducted  by  SRI  International  on  the  generation  and  reflective/absorptive 
properties  of  man-made  atmospheric  pressure  plasma.  The  initial  thrust  was 
to  understand  plasma  generation  in  the  earth's  atmosphere  from  sea  level  to 
300,000  ft.  SRI  thus  investigated  plasma  lifetime;  or  more  exactly  the 
deionization  solution  for  a  tenuous  plasma.  That  investigation  resulted  in  an 
air  chemistry  code  specifically  set  up  for  the  earth's  atmosphere  and  later 
extended  to  noble  gases  [1,2].*  Most  significantly  SRI  derived  an  estimate  of 
plasma  lifetime  as  a  function  of  gas  composition,  pressure,  and  electron  number 
density.  Lifetime  is  important  because  the  power  required  to  sustain  a 
man-made  tenuous  is  inversely  proportional  to  lifetime.  If  the  lifetime  is  long 
enough  the  the  power  required  is  low  and  plasma  generation  becomes  feasible. 

Motivation  for  this  research  was  twofold.  Man-made  plasma  in  the 
earth's  atmosphere  could  be  used  to  (1)  reflect  electromagnetic  waves,  i.e.,  an 
artificial  ionospheric  mirror  [3,4],  or  (2)  absorb  electromagnetic  waves  [5-7]. 

The  cold  collisions!  plasma  model  [8,9]  and  theory  for  reflection  from  an 
Epstein  profile  [10-13]  provided  a  tractable  and  computationally  feasible  model 
for  estimating  electromagnetic  effects.  Estimates  of  the  reflection  coefficient  at 
270,000  ft  [14;  copy  attached]  quantified  the  importance  of  a  sharp  transition 
from  free  space  to  plasma  for  the  ionospheric  mirror  application. 

The  power  reflection  coefficient  for  backscatter  from  a  plasma  was 
computed  and  found  to  be  significantly  lower  than  that  for  commercially 
available  absorber  [15].  The  concept  of  a  cloaking  device,  i.e.,  a  means  to 
reduce  electromagnetic  backscatter,  using  a  collisional  plasma  has  a  relevancy 
to  radar  cross  section  (RCS)  measurements  and  signature  control  [16,17],  A 
plasma  absorber  would  operate  near  the  ambient  atmospheric  pressure  (760 
torr),  be  broadband  from  VHF  through  X-band,  and  have  a  greater  theoretical 
absorption  coefficient  than  any  commercially  produced  carbon-based  absorber. 


*  References  appear  at  the  end  of  this  report. 
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An  important  part  of  this  cloaking  concept  is  estimating  the  power 
required  to  generate  a  plasma  gradient.  The  basic  science  that  affects  the 
practicality  of  a  cloaking  device  includes  (1)  an  efficient  means  of  ionization, 
which  determines  the  energy  required  to  ionize  a  plasma,  md  (2)  the 
deionization  solution  based  on  gas  kinetics,  which  determines  plasma  lifetime 
and  power. 

The  major  thrust  of  this  research  has  been  to  establish  a  model  for 
deionization  in  air  and  noble  gas  plasmas.  Minimization  of  the  power  required 
to  sustain  a  plasma  is  directly  related  to  the  gas  kinetics;  specifically,  the 
momentum-transfer  collision  rate  [18]  and  the  deionization  time  constant  for 
the  specific  gas  mixture  used.  The  air  chemistry  code  developed  as  a  central 
part  of  this  research  quantifies  the  plasma  lifetime  in  air  and  in  noble  gases 
with  trace  impurities  of  air.  The  net  result  has  been  the  identification  of 
helium,  neon,  and  argon  with  impurities  kept  below  100  ppm  as  a  good  (but 
not  minimum)  gas  mixture  for  generating  a  cloaking-device  plasma.  When 
that  research  was  nearly  complete,  attention  was  focused  on  practical  means  of 
ionization . 

Initial  interest  was  on  an  electron-beam  source  with  an  energy  of  100  to 
300  keV,  which  could  generate  the  required  plasma  through  the  mechanism  of 
impact  ionization.  Although  such  a  scheme  works,  has  a  high  input-to-output 
efficiency,  and  has  been  verified  experimentally  by  other  researchers,  it  has 
some  shortcomings.  A  high-energy  electron-beam  source  is  heavy  and  bulky, 
has  a  low  duty  ratio,  and  produces  X-ray  emissions  [19]. 

As  an  alternative,  SRI  has  identified  a  simple,  safe,  compact  means  of 
ionization.  It  is  based  on  ultraviolet  (UV)  photoionization  of  an  organic  vapor 
and  a  confinement  membrane  that  is  both  electrically  transparent  and 
chemically  compatible  with  the  organic  vapor.  This  ionization  method  and 
confinement  membrane  complete  a  system  for  generating  a  cloaking-device 
plasma  that  is  inexpensive  to  construct  and  suitable  for  experimental  validation 
studies. 

Summary  descriptions  of  the  technical  issues  associated  with  air 
chemistry,  broadband  absorption,  and  plasma  generation  are  provided  in  this 
report,  along  with  supporting  documentation. 
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II  TECHNICAL  ISSUES 


This  investigation  has  spanned  nearly  five  years.  Work  on  theory  for 
reflections  and  broadband  absorption  was  completed  last,  year.  The  air 
chemistry  code  was  completed  this  year,  as  was  the  investigation  of  an  organic 
vapor  suitable  for  photoionization  and  a  confinement  membrane. 

A.  Air  Chemistry 

The  predictions  of  plasma  lifetime  in  air  and  in  noble  gases  with 
impurities  build  on  an  existing  data  base  [20-25].  The  range  of  electron 
densities  investigated  was  chosen  to  span  the  applications  of  an  artificial 
ionospheric  mirror  and  a  cloaking  device.  Identification  of  the  limiting 
processes  in  practical  noble  gas  plasmas  [14]  provides  a  theoretical  mode),  for 
prediction  of  attenuation  and  attenuation  bandwidth  as  a  function  of  noble  gas 
species  and  applied  power. 

1.  Additional  Reactions.  The  previously  reported  263  chemical  reactions 
among  43  gas  species  [2]  have  been  increased  to  340  reactions  among  44 
species.  The  addition  of  He(3Si)  provides  a  better  model  for  a  helium  plasma. 
The  additional  two-  and  three-body  reactions  provide  model  He  and  He(3Si) 
reactions  in  plasmas  in  which  their  number  densities  are  appreciable.  The  new 
reactions  (A123-134.  B60-75,  C26-33,  D27-29,  and  F32-62)  are  included  in 
Appendix  A,  which  is  a  complete  tabulation  of  the  reactions  in  the  air 
chemistry  code.  These  reactions  permit  modeling  of  plasmas  in  which  helium 
is  the  dominant  species  and  air  is  an  impurity.  The  air  chemistry  code  can 
(1)  model  the  composition  of  the  earth's  atmosphere  from  sea  level  to  an 
altitude  of  300,000  ft,  or  (2)  model  a  helium  plasma  with  air  as  an  impurity. 

2.  Plasma  Lifetime.  The  air  chemistry  program  has  been  run  to 
determine  the  1/e  lifetime  of  a  plasma  in  the  earth's  atmosphere  as  a  function 
cf  both  electron  number  density  and  alL  ude.  A  plot  of  these  data  is  in 
Reference  14  (copy  attached).  A  complete  listing  appears  in  Appendix  B, 
which  is  a  database  for  validating  the  air  chemistry  code  on  a  new  machine. 

3.  Source  Code.  An  air  chemistry  code  was  developed  on  an  IBM  AT 
computer  equipped  with  a  math  coprocessor  chip.  A  copy  of  the  source  code 
is  attached,  and  the  code  is  briefly  described  in  Appendix  C. 
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B.  Broadband  Absorption 

The  prior  research  on  plasma  absorption  and  backscatter  from  plasma 
gradients  [10—12]  provided  a  solid  theoretical  base.  The  theoretical  model  of  a 
cold  collisional  plasma  [8,9]  is  sufficient  to  predict  electromagnetic  effects  of 
tenuous  atmospheric  plasmas.  The  predictions  of  absorption  and  reflection 
based  on  an  Epstein  profile  provide  an  estimate  of  electrical  characteristics. 

1.  Exact  Modeling.  The  value  of  the  Epstein  profile  gradient  scale  factor 
to  use  for  electron  beam  impact  ionization  or  photoionization  is  approximate. 
The  method  suggested  is  to  plot  the  anticipated  electron  distribution  and  fit 
an  Epstein  profile  to  it.  Clearly,  this  approximate  method  should  be  improved 
by  analytically  or  numerically  computing  electrical  absorption  and  scattering 
for  the  exact  electron  number  density  distribution.  Modeling  the  electrical 
characteristics  is  tractable.  Knowing  the  exact  distribution  is  a  tough 
problem,  because  it  requires  knowing  the  air  chemistry. 

For  example,  a  plasma  generated  in  helium  by  a  photoionization  pulse 
has  a  1/r2  distribution  and  will  deionize  at  a  rate  that  depends  on  the 
electron  concentration.  The  data  in  Appendix  B  for  a  helium  plasma  with  an 
electron  concentration  suitable  for  a  cloaking  device  (103  cm*3  to 
5  x  1011  cm'3)  indicate  that  the  plasma  lifetime  is  a  function  of  electron 
number  density.  Consequently,  the  high-electron-density  portion  of  the  initial 
1/r2  distribution  decreases  in  density  more  quickly  than  the  tail.  And,  the 
Epstein  profile  gradient  scale  factor  flattens  out  (becomes  larger)  as  a  function 
of  time.  Clearly,  the  response  of  a  long  electromagnetic  pulse  will,  in  turn, 
become  time-dependent. 

2.  Bandwidth.  The  Epstein  profile  gradient  scale  factor  [14]  determines 
the  low-frequency  limit,  and  the  electron  momentum-transfer  collision  rate 
determines  the  high-frequency  limit.  The  low-frequency  cutoff  is  inversely 
proportional  to  the  plasma  gradient  scale  factor.  A  plasma  that  varies  from  a 
density  of  n0  to  no/10  in  a  distance  of  2  m  has  a  low-frequency  cutoff  of  150 
MHz.  Lower  or  higher  cutoff  frequencies  can  be  set  by  altering  the  means  of 
ionization  to  produce  the  required  scale  factor. 
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Helium  at  atmospheric  pressure  (14.7  psi)  has  a  momentum-transfer 
collision  rate  of  253  *  109  s'1  [18],  which  produces  a  high-frequency  cutoff  of 
20  GHz.  Argon  and  neon  have  a  momentum-transfer  collision  rate  that  is  an 
order  of  magnitude  lower  and  so  yields  a  2~GHz  cutoff. 

For  a  space-based  cloaking  device  to  have  an  appreciable  high-frequency 
cutoff,  it  must  have  a  confinement  membrane  capable  of  containing  an  internal 
pressure  of  a  few  psi.  An  internal  pressure  of  14.7  psi  of  helium  would  have 
a  high-frequency  cutoff  of  20  GHz,  but  14.7  psi  will  rupture  a  thin  membrane 
several  meters  in  diameter  [26].  By  reducing  the  internal  pressure  by  a  factor 
of  20  to  0.75  psi,  a  high-frequency  bandwidth  of  1  GHz  could  be  realized. 
Consequently,  the  thin-shell  confinement  vessel  described  by  Watters  and 
Vidmar  [26,  copy  attached]  would  be  suitable  for  space  experimentation  at 
frequencies  below  1  GHz. 

C.  Plasma  Generation 

Two  means  of  ionization,  among  numerous  promising  techniques,  have 
been  investigated  in  detail.  High-energy  electron  beam  impact  ionization  and 
photoionization  were  investigated  because  both  can  readily  generate  a  plasma 
with  a  gradient  from  a  single  source.  A  suitable  membrane  material  to 
confine  the  plasma  without  introducing  impurities  or  acting  as  an  electro¬ 
magnetic  scatterer  is  also  necessary.  Significant  results  on  these  aspects  of 
plasma  generation  follow. 

1.  Impact  Ionization.  Electron  impact  ionization  is  a  well-developed 
means  of  generating  a  plasma  in  a  gas.  The  production  rates  for  ionization 
products  and  metastables  in  air  due  to  electron  bombardment  are  known.  A 
production  rate  table  that  approximates  the  generation  of  ionization  by¬ 
products  has  been  incorporated  in  subroutine  XRAY  (in  AIRLFOR). 

Metastable  species  and  other  ionization  by-products  are  important  in  modeling 
conductivity  in  the  late  afterglow  of  an  impulsively  generated  plasma. 

As  a  practical  means  of  powering  a  cloaking  device,  impact  ionization  is 
an  efficient  process.  But,  it  has  some  detracting  features.  First,  the  electron 
must  have  an  energy  over  100  keV  to  penetrate  a  meter  or  so  into  a  gas. 

This  potential  is  sufficient  to  generate  X-rays  [19]  that  constitute  (or  are 
perceived  by  cautious  individuals  to  constitute)  an  unacceptable  health  risk. 
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Second,  a  small,  lightweight,  inexpensive,  nigh- repetition-rate  X-ray  source  is 
not  commercially  available.  A  CW  source  is  more  desirable  but  has  profound 
transmission  window  heating  problems  that  limits  the  time  of  CW  operation. 
Thus,  a  simpler  ionization  technique  was  investigated  as  described  below. 

2.  Photoionization.  UV  ohoton  sources,  4  -  400  nm,  are  commercially 
available  [27],  and  10%  efficiency  at  172  nm  (7.22  eV)  is  possible  [28]. 
Consequently,  a  substance  with  an  ionization  potential  up  to  6  eV  could  be 
ionized  with  a  UV  flash  lamp  using  a  quartz  or  magnesium  fluoride 
transmission  window.  Anderson  [29]  found  that  the  organic  molecule 
tetrakis(dimethylamino)ethylene  (TMAE)  is  oue  of  many  substances  that  can 
be  used  as  a  seed  gas  for  ohotoionizatiou  with  a  UV  flash  lamp. 

Rewick  et  al.  [30]  have  investigated  the  material  compatibility  of  TMAE 
with  common  plastics  and  found  that  Kapton,  Mylar,  and  polyethylene  can  be 
used  to  contain  TMAE  vapor.  The  vapor  pressure  and  cross  section  for 
photoionization  of  TMAE  [29,31]  are  sufficient  for  generating  plasma  with 
dimensions  from  a  few  centimeters  to  many  meters.  Electromagnetic  validation 
measurements  could  be  conducted  using  this  means  of  ionization  if  a  suitable 
confinement  membrane  could  be  constructed. 

3.  Confinement  Membrane.  Two  common  materials,  Kapton  and  Mylar, 
are  candidates  for  the  confinement  membrane  for  a  cloaking  device.  Kapton 
and  Mylar  have  high  tensile  strength  and  a  low  dielectric  constant,  which 
makes  their  backscatter  reflection  coefficients  low.  Watters  and  Vidmar  [26] 
demonstrated  that  2— mil  Mylar  can  be  used  as  a  stressed-skin  conical  structure 
(32-in.  diameter,  96— in.  height,  0.75— psi  relative  pressure)  to  support  targets 
and  antennas  during  electrical  measurements.  This  structure  is  ideal  for 
plasma-absorber  containment  because  of  its  low  backscatter  and  1-m3  volume. 

For  a  cloaking  device,  a  0.25-mil  membrane  of  Kapton  or  Mylar  is 
adequate  for  an  internal  pressure  of  0.1  psi.  Polyethylene  would  have  to  be 
2--mil  thick  and  would  increase  backscatter  from  the  skin  by  18  dB  compared 
with  Mylar.  A  thin  Mylar  membrane  becomes  an  inconsequential  electro¬ 
magnetic  scatterer  below  S-band,  where  80  dB  of  plasma  absorption  is 
predicted  before  reflections  from  the  confinement  membrane  or  plasma  gradient 
become  a  backscatter  issue. 
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Ill  CONCLUSIONS 


After  researching  the  air  chemistry,  electromagnetic  properties,  and 
plasma  generation  mechanisms,  SRI  concludes  that  a  plasma  cloaking  device 
will  work  and  is  feasible.  Additional  basic  research  issues  remain  to  be 
resolved  on  the  cloaking  device  concept  suggested  in  this  report.  Researchers 
in  the  future  could: 

•  Refine  the  air  chemistry  model 

•  Discover  better  seed  gas  mixtures 

•  Develop  more  efficient  ionization  techniques 

•  Quantify  time-  and  source-dependent  electron  density  profiles 

•  Model  scattering  from  time-dependent  density  profiles 

•  Develop  better  confinement  membranes 

•  Discover  ways  to  circumvent  or  slow  rapid  deionization  in  air 

•  Measure  scattering  and  piasma  absorption  coefficients. 

In  addition,  it  is  altogether  possible  that  additional  cloaking  device 
concepts  based  on  different  physical  principles  exist  and  are  simpler. 

The  cloaking  device  concept  advanced  in  this  report  could  be  constructed 
using  helium  or  argon,  seeded  with  TMAE,  photoionized  with  a  UV  flash 
lamp,  and  confined  by  a  Mylar  membrane.  Scattering  and  absorption 
experiments  could  then  be  conducted  using  plane-wave  illumination  with 
free-space  to  plasma  boundary  conditions.  This  confinement  membrane 
provides  a  way  to  generate  a  large  volume  of  plasma  and  is  easy  to  fabricate. 
A  cylindrical  or  conical  confinement  membrane  could  also  be  used  to  generate 
a  plasma  antenna  [32] . 

Applications  of  plasma  absorber  are  those  for  which  carbon-based 
absorber  is  currently  used.  For  applications  in  which  the  reflection  from  the 
tines  of  pyramid  or  wedge  absorber  is  detrimental  or  the  net  absorption  is 
insufficient,  a  plasma  absorber  offers  an  alternative  with  as  much  as  a  40-dB 
improvement.  The  absorption  per  wavelength  of  a  plasma  absorber  is  superior 
to  that  of  a  carbon-based  absorber  and  so  offers  a  compact,  lightweight 
alternative. 
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In  addition  to  affording  excellent  reflection  and  absorption  properties,  a 
plasma  absorber  is  switchable.  This  opens  up  the  possibility  for  its  use  as  a 
high-power  absorber  in  transmission  lines  and  antenna  feeds  to  suppress  ringing 
and  reverberation  in  high-power  pulsed  radars  and  sounders.  The  scattering, 
diffraction,  and  radiation  lobe  characteristics  of  antennas  could  be  altered  by 
appropriate  application  of  plasma  absorber.  These  and  other  potential 
applications,  such  as  an  ionospheric  mirror  and  plasma  antenna,  provided  the 
basic  motivation  for  this  research  contract. 
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Appendix  A 

ATMOSPHERIC  REACTION  RATES 


The  numeration  system  for  the  reactions  in  this  Appendix  continues  the 
system  adopted  in  previous  SRI  documentation  for  this  project  [1,2].  The 
reactions  have  been  reorganized  to  group  sequentially  all  the  reactions  relating 
to  a  single  species.  The  index  number  preceding  a  reaction  (such  as  All  or 
F22)  corresponds  to  the  vari  .ble  name  All  and  F22  ir  the  air  chemistry  code 
Because  the  reactions  were  added  sequentially  withou  ogard  for  species  and 
then  reorganized,  the  numeration  in  the  following  lists  is  not  sequential.  The 
340  reactions  below  relate  to  the  44  species  in  Table  A-l.  All  reaction  rates 
were  extracted  from  Bortner  and  Bauer  [21]  except  for  the  ones  followed  by 
BA,  which  were  extracted  from  Baulch  et  al.  [33].  Reactions  that  have  not 
yet  been  measured  have  been  assigned  a  typical  value  extracted  from  Bortner 
and  Bauer  [21].  The  units  for  two-body  and  three-body  reactions  are  cm3/s 
and  cm6/s,  respectively. 


Table  A-l 

SPECIES  IN  SIMULATION  CODE 


Neutral  and  Metastable 

Positive 

Negative 

H,  He,  N,  0 

n2,  o2,  no,  oh 

°3 

C02,  HjO,  N20,  N02 
He(3Sj) 

0(‘S),  O('d) 

N(2P°),  N(2D°) 

02fa‘Ag),  02(b1E+) 

N+,  N+,  N+ 

o+,  o+  0+ 

He+,  He2 

h2o+,  h3o+ 

NO+ 

0+.H20 

h3o+»h2o 

h3o+*oh 

e 

0 , 02,  Og,  04 

o2.h2o 
co~  co; 

co“»h2o,  co~*h2o 

A.  NEGATIVE-ION  REACTIONS 


Reaction 


Al. 

e 

4 

°2 

-4 

0  + 

0 

A2. 

e 

+ 

°3 

0  4- 

°2 

A3. 

e 

4 

o2 

4 

t 

o 

4 

°2 

A4. 

e 

4 

°2 

4 

n2- 

°2 

4 

N2 

A5. 

e 

4 

°2 

4 

h20 

-  0 

2 

+  H2° 

A6. 

e 

4 

°2 

4 

co2 

-  0 

2 

+  C02 

A7. 

e 

+ 

°2 

4 

He 

°2 

4 

He 

A132.  e  4  He(3Sj)  ^  He  4  e 


__ _ _ Rate _ 

1.0E(-16) 

9.0E(-12) 

1.4E(-29)  (300/Te)  exp(-600/Te) 
1.0E(-31) 

1.4E(— 29) 

3.3E(-30) 

1.0E(-31) 

7.0E(-10)(Te/300)1/2 


A8. 

O' 

'  -f 

0  -*  02  +  e 

2.0E(— 10) 

A9. 

O' 

■  4 

N  -»  NO  4  e 

2.0E(-10) 

A10. 

O' 

■  4 

°2  “*  °3  +  e 

5.0E(-15) 

All. 

0” 

'  4 

02(a  Ag)  -*  Og 

4 

e 

3.0E(-10) 

A12. 

O' 

'  4 

02(»lAg)  -  0- 

4 

0 

l.OE(-lO) 

A13. 

O' 

4 

°3  °3  4  0 

5.3E(-10) 

A14. 

O' 

4 

-♦  N20  4  e 

1.0E(-14) 

A15. 

0“ 

+ 

NO  -*  N02  4  e 

2.5E(-10)  (300/Te)0-8 

A16. 

0“ 

4 

^2  +  ^2  ^3 

4 

°2 

1.0E(-30) 

A17. 

O' 

4 

02  +  C02  -»  COg 

+  Or, 

3.1E(-28)  (300/Te) 

A18. 

0“ 

4 

02  4  He  0^ 

4 

He 

1.0E(-30) 

A  64. 

o' 

4 

02(b1E^)  -  0 

4 

°2  +  e 

l.OE(-lO) 

A76. 

0“ 

4 

O^D)  -*  0  4 

0 

4-  e 

l.OE(-lO) 

A89. 

o' 

4 

N(2D°)  -*  0  4 

N 

4  e 

I.OE(-IO) 
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A105. 

0~  4-  N(2P°)  -.0  4-  N  4-  e 

I.OE(-IO) 

A122. 

0“  4-  0(1S)  -  0  +  0  4-  e 

l.OE(--lO) 

A123. 

0“  4-  H e(3Sj)  -  0  +  He  +  e 

3.0E(-10) 

A19. 

O2  4-  0  -»  0  4-  Og 

1.5E( -10) 

A20. 

0“  4-  N  -  N02  +  e 

3.0E(-10) 

A21. 

0~2  +  N  ->  0~  +  NO 

LOE(--IO) 

A22. 

02  4-  02  -»  02  4-  02  4-  e 

2.7E(-10)  (Tn_/300)0-5 
U2 

x  exp(-5590/To_) 

A23. 

02  4-  02(a  Ag)  -*  02  4-  00  4-  e 

2.0E(-19) 

A124. 

02  4-  He(  Sj)  -*  02  4-  He  4"  e 

3.0E(-10) 

A24. 

°2  4-  °3  °2  +  °3 

4.0E(-10) 

A25. 

°2  +  N2  °2  +  N2  +  e 

1.9E(-12)  (Tn_/300)L5 
u2 

*  exp(-4990/Tn_) 
u2 

A2(>. 

°2  4-  °2  4-  °2  -  °;  +  02 

3.5E(— 31)  (300/Tair) 

A27. 

0~  +  02  4-  H20  -  02*H20  +  02 

3.0E(--28)  (300/Tair) 

A28. 

02  4 -  02  4-  COr,  -♦  C04  4-  02 

4-7E(— 29)  (300/Taif) 

A29. 

02  4-  02  4-  He  -*  0^  4-  He 

3.5E(-31)  (300/Tair) 

A63. 

0~  4-  02(b1E^)  -  02  4-  02  4-  e 

l.OE(-lO) 

A75. 

0~  +  0(1D)  -  02  +  0  +  e 

l.OE(-lO) 

A88. 

0“  4-  N(2D°)  -»  02  4-  N  +  e 

l.OE(-lO) 

A104. 

0“  4-  N(2P°)  -  02  4  N  4-  e 

l.OE(-lO) 

A 120. 

0“  +  0(!S)  -  0“  4  0  4-0 

l.OE(-lO) 

A121. 

02  4-  0(1S)  -  02  4-  0  4-  e 

l.OE(-lO) 
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A30. 

1  ec 
C 

4 

0  -■*  02  4 

°2  +  ' 

3 

LOE(-ll) 

A31. 

°3 

4 

0  -  0~  4 

°2 

3.2E(-10) 

A32. 

°i 

4 

co2  -•»  co~ 

4  02 

5.5E(~10)  (300/T 

A54. 

°5 

4 

02(b>E+)  . 

*  0“  4 

°2  + 

°2 

l.OE(-lO) 

A65. 

°5 

4 

0(1D)  -4  O' 

4  02 

4  0 

l.OE(-lO) 

A77. 

°3 

+ 

N(2D°)  -4  0 

“40 

2  +  N 

I.OE(-IO) 

A78. 

°3 

4 

N(2D°)  -  0 

3  +  N 

4  e 

l.CE(-lO) 

A90. 

°3 

4 

N(2P°)  -  0 

“  4  o; 

?  +  N 

l.OE(-lO) 

A9i. 

+ 

N(2P°)  -4  0 

9  4  0 

4  N 

I.OE(-IO) 

A92. 

+ 

N(2P°)  -4  0 

3  +  N 

4  e 

1.0E(~10) 

A106. 

°S 

4 

O^S)  -  0~ 

4  Or, 

4  0 

l.OE(-lO) 

A107. 

°3 

+ 

0(iS)  -  0” 

4  0 

+  0 

I.OE(-IO) 

A108. 

°; 

4 

0(1S)  -  03 

4  0  4  e 

I.OE(-IO) 

A 125. 

°3 

4 

He^Sj)  0 

'2  4  0 

4  He 

4  e 

3.0E(~10) 

A33. 

4- 

0  -»  Og  4  02 

4.0E(-10) 

A34. 

°i 

4 

°2  ■*  °2  + 

°2  + 

°2 

2.7E(-14) 

A35. 

°4 

4 

°3  “*  °3  + 

°2  + 

°2 

3.0E(~10) 

A36. 

°4~ 

4 

co2  -*  CO^ 

4  02 

4.3E(~10) 

A37. 

°I 

4 

H20  -4  0~«H23  4 

°2 

1.4E(-9) 

A49. 

4 

02(a1Ag)  - 

°2  + 

°2  + 

°2 

l.OEf-10) 

A59. 

4 

02(b1E+)  - 

°2  + 

°2  + 

°2 

l.OE(-lO) 

A60. 

°i 

4 

o2(blz+)  - 

°2  + 

°2  + 

°2  4  e 

l.OE(-lO) 

A71. 

°; 

4 

O^D)  -»  0~ 

4  02 

4  0 

l.OE(-lO) 

A?2. 

4 

0(1D)  00 

4  O9 

4  0 

4  e 

3.0E(-1C) 

0.49 
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A84. 

Q~  4  N(2D°)  -  OJ  4  02  4-  N 

l.OE(-lO) 

A85. 

0~  4  N(2D°)  -  02  +  02  +  N  4  e 

l.OE(-lO) 

A99. 

Q~  4  N(2P°)  -*  02  +  02  4  N 

l.OE(-lO) 

AlOO. 

0"  4  N(2?0)  -  02  4  02  4  N  4  e 

l.OE(-lO) 

A115. 

0"  4  0(!S)  -  o2  +  o2  +  o 

l.QE(-lO) 

A116. 

0^  4  0(*S)  -*  02  4  02  +  0  4  e 

1.0E(~10) 

A 126. 

+  He(3S1)  -»  ©2  4-  ©2  4  He  4  e 

3.0E(-10) 

A33. 

^2*^2^  ^3  ”*  ^3  "b  4  H20 

2.3E(-10) 

A39. 

o2*h2o  4  co2  -*  co~  4  h2o 

5.8E(-10) 

A51. 

0”«H20  4  02(a1Ag)  -  0~  4  K20  4  02 

l.OE(-lO) 

A55. 

0~«H20  4  02(blEj)  -  0~  4  H20  4  02 

1.0E^-10) 

A  56. 

0“«H20  4  02(b1Sg)  -  02  4  H20  4  02  4  e 

l.OE(-lO) 

A66. 

°2*H2°  +  °(1d)  ^  °2  +  H20  4  0 

1.0E(  -10) 

A67. 

02«H20  4  O^D)  -*  02  4  H20  4  0  4  e 

l.OE(-lO) 

A79. 

C2«H20  4  N(2D°)  -  0~  4  H20  4  N 

l.OE(-lO) 

A80. 

02«H00  4  N(2D°)  -  02  4  H20  4  N  4  e 

l.OE(-lO) 

A93. 

0~*H20  4  N(2P°)  -  02”  4  H20  4  N 

l.OE(-lO) 

A94. 

0~»K20  4  N(2P°)  -  02  4  H20  4  N  4  e 

l.OE(-lO) 

A109. 

02*H20  4  O^S)  -  0~  4  Ii20  4  0 

LOE(-IO) 

AllO. 

02*H20  t  0(lS)  -»  02  4  H20  4  0  4  e 

l.OE(-lO) 

A127. 

02*H20  4  He(3Sj)  -  02  4  H20  4  He  4  e 

3.0E(-10) 

A40. 

CO3  4  0  -*  02  4  C02 

I.IE(-IO) 

A41. 

CO3  4  02  4  HgO  C02*H29  4  02 

1.0E(~28)  (300/Tair) 

A133. 

CO"  4  He  4  H20  -  C03*H20  4  He 

1.0E(-28)  (300/Taif) 
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A68. 

co- 

4 

0(2D)  -  0"  -1-  C02  4  0 

I.OE(-IO) 

A81. 

co- 

4 

N(2D°)  0"  4  C02  4  N 

I.OE(-IO) 

A95. 

co: 

4 

N(2P°)  -  0"  4  C02  4  N 

l.OE(-lO) 

A96. 

C°3 

4 

N(2P°)  4  0  4  C02  4  N  ■ 

f  e 

l.OE(-lO) 

A  .1 1 1 . 

C°3 

+ 

0(1S)  4  0"  4  C02  4  0 

l.OE(-lO) 

A112. 

CO3 

4 

0(1S)  -  0  4  C02  4  0  4 

e 

l.OE(-lO) 

A128. 

CO¬ 

4- 

He(3S1)  4  C02  4  0  4  He 

4  e 

3.0E(-1C) 

A42. 

CO" 

4 

0  4  CO3  4  02 

1.5E(— 10) 

A43. 

co- 

4 

°2  -  °4  4  C02 

2.0E(— 14) 

A44. 

co¬ 

+ 

O3  4  O3  4  C02  4  02 

1.3E(-10l 

A45. 

co, 

4 

4  H20  4  C04*H20  4 

N2 

5.0E(-29)  (300/Tajr) 

A 134. 

co¬ 

4 

He  4  H20  4  C0"*H2O  4 

He 

5.0E(— 29)  (300/T^) 

A46. 

co" 

4 

H20  4  02»H20  4  C02 

2.0E(-10) 

A50. 

CO" 

4 

02(a1Ag)  4  02  4  C02  4 

°2 

I.OE(-IO) 

A57. 

C°4 

4 

°2(blsJ)  4  0~  4  C02  4 

°2 

l.OE(-lO) 

A58. 

CO" 

4 

O^E*)  4  02  4  C02  4 

02  +  e 

l.OE(-lO) 

A69. 

CO" 

4 

O^D)  4  0"  4  C02  4  0 

l.OE(-lO) 

A70. 

COJ 

4 

o^d)  -  o2  4  co2  4  0 

4  e 

l.OE(-lO) 

A82. 

co¬ 

4 

N(2D°)  4  O'  4  C02  4  N 

l.OE(-lO) 

A83. 

co" 

4 

N(2D°)  4  02  4  C02  4  N 

4  e 

l.OE(-lO) 

A97. 

CO" 

+ 

N(2P°)  4  0“  4  C02  4  N 

l.OE(-lO) 

A98. 

co4 

4 

N(2P°)  4  02  4  C02  4  N 

4  e 

I.OE(-IO) 

A113. 

co- 

4 

0(!S)  4  0“  4  C02  4  0 

l.OE(-lO) 

All  4. 

co¬ 

4* 

0(XS)  4  02  4  C02  4  0 

4  e 

l.OE(-lO) 

A129. 

co" 

He(3Sj)  4  C02  4  02  4  He  4  e 

3.0E(-10) 
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A47, 

C03»H20 

+ 

n2  ^  CO"  +  h20  +  N2 

3.0E(-14) 

A52. 

CCyHgO 

+ 

°2(a1Ag)  -  CO"  +  H20  +  02 

l.OE(-lO) 

A61. 

C03»H20 

+ 

0^+)  -  CO3  +  h20  +  02 

l.OE(-lO) 

A73. 

C03*Ho0 

+ 

O^D)  CO"  +  H20  +  0 

LOE(-IO) 

A86. 

co“*h2o 

+ 

N(2D°)  -*  CO~  +  H20  +  N 

I.OE(-IO) 

AlOl. 

C0g*H20 

+ 

N(2P°)  -4  OO3  +  H20  +  N 

l.OE(-lO) 

A102. 

COg-I^O 

4- 

N(2P°)  -4  C02  +  H20  +  N 

l.OE(-lO) 

4-  0  +  e 

All  7. 

coj»h2o 

+ 

0(1S)  -4  CO'  +  h2o  4-  0 

I.OE(-IO) 

A118. 

C0g*H20 

+ 

0(1S)  -4  co2  +  H20  +  0 

1.0E(~10) 

4-  0  +  e 

A130. 

C03«H20 

+ 

He(3S1)  -  C02  4-  0  +  H20 

3.0E(-10) 

4-  He  +  e 

A48. 

co~«h2o 

+ 

n2  -*  co~  +  h2o  +  n2 

1.0E(— 14) 

A53. 

co~*h2o 

+ 

02(a1Ag)  -4  CO~  4-  H20  4-  02 

l.OE(-lO) 

A62. 

C0~*H20 

+ 

°2(blSg)  -4  CO"  4-  H20  4-  02 

I.OE(-IO) 

A74. 

C0"*H20 

+ 

0(1D)  -*  CO"  4-  H20  +  0 

l.OE(-lO) 

A87. 

co^»h2o 

+ 

N(2D°)  -4  CO~  +  H20  +  N 

I.OE(-IO) 

A103. 

co~»h2o 

+ 

N(V)  -4  CO~  f  H20  +  N 

l.OE(-lO) 

A119. 

co~»h2o 

+ 

o^s)  -4  CO"  +  h2o  +  0 

I.OE(-IO) 

A131. 

co~»h2o 

+ 

He(3S1)  -  C02  +  02  +  H20 

3.0E(-10) 

4-  He  +  e 
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B,  POSITIVE-ION  REACTIONS 

Reaction 


Rate 


Bl.  N4 
B2.  N4 
B3.  N4 

B4.  N+ 
B5.  N4 
B6.  N4 
B64.  N4 

B7.  N4 
B33.  N4 
B37.  N4 
B42.  n£ 
B47.  N^ 
B52.  N4 
B57. 

B63.  N4 

B8.  O4 
B9.  O4 
BIO.  O4 
B73.  O4 


+  O2  ■*  N  +  04 

3.0E(-10) 

O2  -*  NO4  ■+•  0 

2.8E(— 10) 

+  H20  -»  H204  +  N 

2.6E(-9) 

+  0  -  NO4  +  N 

1.3E(— 10)  (300/T^)0  46 

+  N2  +  N2  -  N+  +  N2 

5.0E(-29) 

+  N2  +  He  N4  +  He 

5.0E(-29) 

-f  He^Sj)  -*  N4  +  N  +  He 

I.OE(-IO) 

+  °2  °2  4  N2  4  N2 

4.0E(-10) 

+  -  N+  +  N2  +  02 

l.OE(-lO) 

+  02(b1E4)  -  N4  +  N2  +  02 

1.0E(-10) 

+  O^D)  -  N4  +  N?  +  0 

l.OE(-lO) 

+  N(2D°)  ->  n4  +  n2  +  n 

1.0E(-10) 

+  N(2P°)  -  N4  +  N2  +  N 

l.OE(-lO) 

+  O^S)  -  N4  +  N2  +  0 

1.0E(-10) 

+  He(3Sj)  -*  N+  +  N  +  N2  +  He 

1.0E(-10) 

+  o2  04  +  0 

2.0E(-11)  (300/T^p0-4 

+  H20  -»  H204  +  0 

2.33E(-9) 

4  N2  4  N2  N0+  +  N  4-  N2 

6.0E(-29)  (300/T^j.)2 

+  N2  +  He  -*  NO4  +  N  +  He 

6.0E(-29)  (300/T  •  )2 

dl  L 
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Bll. 

+ 

N  -*  N0  + 

+  o 

l..2E(-10) 

B12. 

°2 

+ 

02  +  02 

"*  O4  + 

°2 

2.8E(-30) 

B13. 

°2 

+ 

02  +  He 

-  0+  + 

He 

2.8E(~30) 

B14. 

°2 

+ 

H20  +  N, 

2  -  Oj.H20 

+ 

N2 

2.8E(-28)  (300/Tair)i 

B62. 

°2 

+ 

He(3Sj)  -4 

0+  +  0 

'  + 

He 

l.OE(-lO) 

B74. 

+ 

H20  +  He  -  0^«H20 

+ 

He 

2.8E(-28)  (300/Tajr)5 

B15. 

°4 

+ 

O  "  °2  +  °3 

3.0E(-10) 

B16. 

°4 

+ 

O^a'Ag) 

-0+  + 

°2 

+ 

°2 

LOE(-IO) 

B17. 

°4 

+ 

N2  "*  °2 

+  °2  + 

N2 

2.6E(-5)  (300/Tajr)4'2 

x  exp(~5400/Tajr) 

B18. 

°4 

+ 

H20  4  0^«H20  + 

°2 

L5E(-9) 

B39. 

°4 

+ 

-0+  + 

°2 

+ 

°2 

LOE(-IO) 

B44. 

°4 

+ 

0(XD)  -  0+  +  02 

4*  ^ 

0 

I.OE(-IO) 

B49. 

°4 

+ 

N(2D°)  - 

0+  +  Of 

)  + 

N 

I.OE(-IO) 

B54. 

°4 

+ 

N(2P°)  - 

°2  +  °5 

>  + 

4 

N 

LOE(-IO) 

B59. 

°4 

+ 

0(1S)  -  0 

2  +  °2 

+  0 

l.OE(-lO) 

B61. 

°4 

He(3Sj)  -* 

0+  +  0 

+ 

°2 

+  He 

l.OE(-lO) 

B19. 

He+ 

+ 

n2-n+ 

+  N  +  He 

6.0E(-1C) 

B20. 

He+ 

4- 

n2-n+ 

+  He 

e.oE(-io) 

B21. 

He+ 

+ 

o2-o+ 

+  0  +  : 

He 

6.0E(-10) 

B22. 

He+ 

+ 

°2  "  °2 

+  He 

6.0E(— 10) 

B23. 

He+ 

He  +  He 

“  He2  + 

He 

1 . 1E(— 31) 
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B24. 

He2  +  N2  ■*  N2  +  He  +  He 

1.2E(-9) 

B60. 

He 2  4  02  -  °2  4  He  4  He 

1.2E(-9) 

B65. 

He 2  4  He(3Sj)  -*  He"1"  4-  He  +  He 

I.OE(-IO) 

B68. 

N0+  +  He(3S5)  -  N+  4  0  4  He 

5.0E(-11) 

B71. 

N0+  4  He(^Sj)  -*  04  +  N  +  He 

5.0E(-11) 

B25. 

o+*h2q  4  o2  -+  0+  4-  H90 

9.4E(-14) 

B26. 

°2*H2°  +  °2(alAg)  -  °2  +  H2°  +  °2 

l.OE(-lO) 

B27. 

°2*H2°  4'  N0  N0+  +  h2°  +  °2 

I.OE(-IO) 

B28. 

o2«»h2o  +  h2o  h3o+«oh  +  o0 

1.0E(-9) 

B38. 

02.H20  4-  02(b  ^g)  "*  ^2^  4  4  02 

I.OE(-IO) 

B43. 

0+*H20  4  O^D)  -oj  4  H20  4  0 

l.OE(-lO) 

B48. 

°2*H2°  +  N(2D°)  -*  0^  4-  H20  4-  N 

l.OE(-lO) 

B53. 

0+«H20  4-  N(2P°)  -*  02  4-  H20  f  N 

l.OE(-lO) 

B58. 

04  *H20  4-  0(h)  0‘J  4  h2o  4  0 

l.OE(-lO) 

B69 

°2*H2°  +  He(3Sj)  -  0+  4  0  4  H20  4  He 

l.OE(-lO) 

B99. 

H20+  4  02  -*  0*2  4  H20 

2.0E(-10) 

B30. 

H20+  4  H20  -  H30+  4  OH 

1.8E(-9) 

B66. 

H204  4  He(3S1)  -  0+  4  2H  4  He 

I.OE(-IO) 

B31. 

H30+  4  H20  4  N2  m  H30+«H90  4  N2 

3.4E(-27) 

B67. 

H30+  4  He(3S1)  ->  0+  4  3H  4  He 

l.OE(-lO) 

B75. 

H30+  4  H20  4  He  -»  H30+*H20  4  He 

3.4E(~27)  (300/Tair) 
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B32.  HgO^OH  4  H20  -  HgO+^G  4  OH  L0E(-9) 

B34.  H3Of*OH  4  02(a1Ag)  -  H30+  4  OH  +  02  I.OE(-IO) 

B36.  H30+*0H  4  02(h1Ej)  -  H30+  4-  OH  +  02  I.OE(-IO) 

B41 .  H30+-0H  4-  0(3D)  -*  H304  4  OH  4-  0  l.OE(-lO) 

B46.  H30+«0H  4  N(2D°)  ->  H30+  +  OH  4  N  l.OE(-lO) 

B51.  H30+*0H  4  N(2P°)  -  H30+  4  0H4N  l.OE(-lO) 

B56.  H30+«0H  +  0(3S)  -  H30+  4  OH  4  O  l.OE(-lO) 

B72.  H304«0H  4  He(3S1)  -  0+  4  3H  4  OH  4  He  I.OE(-IO) 

B35.  H304'*H20  4-  02(b1S4)  -  H3Of  4-  H?0  4  02  l.OE(-lC) 

B40.  H30+«H20  4-  0(3D)  H30+  4  H20  4-  O  LOE(-IO) 

B45.  H30f«H20  4  N(2D°)  -  HgO4  4-  HgO  +  N  l.OE(-lO) 

B50.  H30+.H20  4-  N(2P°)  ->  H304  4-  H20  4-  N  l.CE(-lO) 

B55.  H30  f*H20  4-  0(3S)  HgO4  4-  HgO  4-  O  LOE(-IO) 


B70.  HjO+^O  4  He(3Sj)  -  0+  j-  3H  4-  H20  4-  He  LOE(-IO) 
C.  NEUTRAL-SPECIES  REACTIONS 


Reaction  _  Rate 


Cl.  N  4-  N  4-  N2  -  N2  4-  N2 

7.6E(-34)  exp(500/Tair) 

C31.  N  4  N  4  He  ■*  N2  4  He 

7.6E(-34)  exp(500/Tajr) 

C2.  0  4  0  4  N2  -  02  4  N2 

3.0E(-33)  (300/T  .  )2’9 

air7 

C32.  0  4-  G  4"  He  -*  02  4-  He 

3.0E(-33)  (300/T  .  )2  9 

ail 

C3.  0  4-  02  4-  N2  -  03  4-  N2 

5.5E(-34)  (300/T  .  )2’6 
air 

C33.  0  4-  02  +  He  ^  03  +  He 

5.5E(-34)  (300/T^p2'6 
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C7. 

O3  4-  02(b1Eg)  ^  02  +  02  +  0 

2.2E(-11) 

BA 

C8. 

O3  +  0(1D)  -  02  +  02 

1.2E(— 10) 

BA 

C9. 

O3  f-  0(!D)  -4  o2  +  0  +  0 

1.2E(— 10) 

BA 

CIO. 

O3  +  N(2D°)  -4  02  +  N  +  0 

l.OE(-lO) 

Cll. 

O3  +  N(2P°)  4  02  +  N  +  0 

l.OE(-lO) 

C12. 

O3  +  O^S)  -  02  +  0  +  0 

l.OE(-lO) 

C13. 

N20  +  0(1D)  -  N2  +  0  +  0 

l.OE(-lO) 

C14. 

N20  +  O^D)  -4  2NO 

6.6E(-11) 

BA 

C15. 

n2o  +  N(2D°)  4  n2  +  N  +  0 

l.OE(-lO) 

C16. 

N20  +  N(2P°)  4  n2  +  N  +  0 

l.OE(-lO) 

Cl  7. 

N20  +  0(XS)  4  N2  +  0  +  0 

l.OE(-lO) 

C18. 

N09  +  N(2P°)  NO  +  0  +  N 

LOE(-IO) 

C19. 

N02  +  0(1S)  -•  NO  +  0  +  0 

l.OE(-lO) 

C4. 

°2(a*Ag)  +  02  -  02  +  02 

2.4E(-18) 

C5. 

02(a1Ag)  +  N2  -  02  +  N2 

1 . 1E(— 19) 

C6. 

02(a1Ag)  +  H20  -  02  +  H20 

1 .5E(— 17) 

C20. 

02(b1E+)  +  N2  4  02  +  N2 

2.0E(-15) 

BA 

C21. 

C2(b  Eg)  +  02  -4  09  +  02 

4.0E(~17) 

BA 

C22. 

+  H20  -  02  +  H20 

4.0E(-12) 

BA 
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C23.  0(1D)  +  N2  -  0  4-  N2 

l.SE(-ll)  exp(107/Tair)  BA 

C24.  CK1!))  +  02  -  0  +  02 

3.2E(-11)  exp(67/Tair)  BA 

C25.  O^D)  +  H20  -  20H 

2.3E(-10)  BA 

C26.  He(3Sj)  +  He(3S1)  -  He+  4  e 

2.0E(-9) 

C27.  He(3S1)  +  He(3S1)  -*  He+  +  He  4-  e 

2.0E(-9) 

C28.  He^Sj)  4-  N2  -  N+  +  N  +  He  4-  e 

l.OE(-lO) 

C29.  He(3Sl;  +  02  -*  0+  +  0  +  He  4-  e 

l.OE(-lO) 

C30.  He^Sj)  4-  H20  -»  H20+  +  He  4-  e 

l.OE(-lO) 

D.  POSITIVE-ION  ELECTRON  RECOMBINATION 

Reaction 

Rate 

Dl.  N+  +  e  +  e-*N  +  e 

7.0E(-20)  (300/Te)4,5 

D2.  N4"  +  e  +  N2  ->  N  -f  N2 

6.0E(-27)  (300/Te)2'5 

1)28.  N+  +  e  +  He  -*  N  +  He 

6.0E(-27)  (300/T  )2,5 

D3.  +  e  +  e  -*  N2  +  e 

7.0E(-20)  (300/Te)4,5 

D4.  N2  4  e  -*  N  -f  N 

2.7E(-7)  (300/Te)0'2 

D24.  4  e  +  e  -*  N2  +  N2  +  e 

7.0E(-20)  (300/Te)4,5 

D5.  N+  4  e  -*  N2  4  N2 

2.0E(-6) 

D6.  0+  4e4e-*04e 

7.0E(-20)  (300/Te)4'5 

D7.  O'  4  e  4  02  -*0  4  02 

6.0E(-27)  (300/Te)2'5 

D29.  0+  4  t  4  He  0  4  He 

6.0E(-27)  (300/Te)2"5 
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D8. 

0+  4-  e  +  e  -*  O2  +  e 

7.0E(-20)  (300/Te)4-' 

D9. 

0+  4-  e  -*  0  4-  0 

2.1E(-7)  (300/TJ0-7 

D25. 

0^  4-  e  -r  e  -*  O2  4*  O2  4-  g 

7.0E(-20)  (300/Te)4  £ 

DIO. 

0+  +  6  -•  0  +  0  +  (>2 

2.0E(  -6) 

Dll. 

He+  +  e  +  e-»He  +  e 

1.0E(~19) 

D27. 

He+  4-  e  +  e  -»  He^Sj)  4-  e 

3.0E(-20)  (300/Te)4'5 

D23. 

He^  +  e  +  e-*He  +  He  +  e 

7.0E(-20)  (300/Te)4  5 

D12. 

He2  +  e  -»  He  4-  He 

L0E(-8) 

D13. 

He^  4-  e  4-  He  ->  3He 

2.0E(~27) 

D14. 

N0+  4-  e  4-  e  -»  NO  4-  e 

7.0E(-20)  (300/Te)4,5 

D15. 

N0+  4-  e  N  4  0 

4.0E(-7)  (300/T  )0'4 

V* 

D16. 

HjO  4-  e  4-  e  -»  H2O  4'  H  4-e 

7.0E(-20)  (300/Te)4-5 

D17. 

H30+  4-  e  -*  H20  4-  H 

1.3E(-6)  (300/Te)10 

D26. 

O^HgO  4  e  4  e  ■*  02  4-  H20  4*  e 

7.0E(-20)  (300/Te)4  5 

D18. 

0^.H20  4-  e  -*  02  4-  H20 

1.5E(-6)  (300/Te)0  2 

D19. 

H30+»0H  4-  e  4-  e  -  H<20  +  OH 

4-  H  4-  e 

7.0E(-20)  (300/Te)4’5 

D20. 

H30+«0H  4-  e  -  H90  +  OH  4-  H 

2.0E(-6)  (300/Te)0-2 

22 


D21.  H3CT»H20  4-  e  +  e  -  2H20  +  H  +  e  7.0E(~20)  (300/TJ45 
D22.  H30+*H20  +  e  -  2H20  +  H  2.8E(-6)  (300/Te)0,15 

E.  TWO-BODY  POSITIVE-ION  NEGATIVE-ION  RECOMBINATION 


Reaction  Rate 


El. 

N+ 

+  0  -  N  f  0 

2.6E(-7) 

(300/Tajr)0'5 

E2. 

N+ 
1N  2 

+  o2  -»  n2  +  o2 

1.6E(-7) 

(30°/Tair)0'5 

E3. 

04' 

+  0"  -  0  +  0 

2.7E(— 7) 

(300/Tajr»0'^ 

E4. 

°2 

+  o  -*  o2  +  0 

1.0E(— 7) 

(3°0/Tair)0'5 

Ej. 

°2 

+  °2  "*  °2  +  °2 

4.2E(-7) 

(SM/Tair)0'5 

E6. 

NO+  +  0  -*  NO  +  0 

4.9E(— 7) 

pOO/T^r)0'5 

E7. 

N0+  -1-  02  -*  NO  +  02 

6.0E(-7) 

(300/T^r)0-5 
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F.  THREE-BODY  CLUSTER-ION  RECOMBINATION 


The  estimate  cited  by  Bortner  and  Bauer  [21]  for  the  three-body 
cluster-ion  recombination  rate  is  l.GE(-25)  (300/T  .  )5/"2  where  the  third  body 

ctll 

is  N«2  or  the  dominant  gas  species.  Because  little  information  is  available  for 
individual  reactions,  the  reaction  rate  for  all  the  following  reactions  is  assumed 
to  be  the  same. 


Reaction 


FI. 

°2 

+ 

®3  +  N2 

202 

4 

0  4 

N2 

F2. 

°2 

4 

4  N2  - 

302 

4 

N2 

F3. 

°2 

+ 

C03  +  N2 

°2 

4 

co2 

4  0 

4* 

n2 

F4. 

°2 

4* 

CO"  4  N, 

*  ^ 

-•  20 

2  4  CO, 

2  +  N2 

F5. 

°2 

4 

°2®H2°  + 

N2- 

202  4 

h2° 

4 

N2 

F32. 

°2 

4 

Oj  4  He  -* 

202 

4 

0  4 

He 

F33. 

°2 

■t* 

O4  4  He  - 

302 

4 

He 

F34. 

°2 

4 

CO"  4  He 

-°2 

4 

co2 

4  0 

4 

He 

F35. 

0+ 

4 

CO"  4  He 

-*  202  4  COf 

>  4  Ke 

m 

F36. 

°2 

4 

o-.h2°  + 

He  - 

202  + 

h2° 

4 

He 

F6. 

°4 

4 

0"  4  N?  - 

20o 

A* 

4 

0  4 

N2 

F7. 

0+ 

u4 

4 

°2  +  N2  “* 

302 

4 

n2 

F8. 

<>4 

4 

O3  4  N0  -» 

202 

+ 

°3  + 

*2 

F9. 

°4 

4 

0“  4  N2  - 

400 

Am 

4- 

W2 
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F37.  0+  4-0  4  He  -»  2O2  4  0  4  He 

F38.  0+  +  O”  4  He  -  302  4-  He 

F39.  0+  4  O3  4  He  -»  202  4  O3  4  He 

F40.  OJiOj+He.  402  +  He 

F10.  N0+  4  O3  4  N2  -  NO  4  02  4  O  4 

Fll.  N0+  +  0“  4  N2  -  NO  4  202  4  N2 

F12.  NO^  4  C0“  4  N2  h  NO  4  C02  4  0  4  N2 

F13.  N0+  4  C0~  4  N2  h  NO  4  C02  4  02  4  N0 

F14.  N0+  4  0~*H20  4  N2  -*  NO  4  HgO  4  O  4  N2 

F15.  N0+  4  C0“*H20  4  N2  -4  NO  4  H20  4  C0<2  4  O  4  N<2 

F16.  N0+  4  CO^HgO  4  N2  -  NO  4  HgO  4  C02  4  02  4  N2 

F41.  N0+  4  O3  4  He  -*  NO  4  Oj  4  0  4  He 

F42.  N0+  4  0”  4  He  -+  NO  4  202  4  He 

F43.  N0+  4  CO3  4  He  NO  4  C02  4  0  4  He 

F44.  N0+  4  C0~  4  He  -  NO  4  C02  4  02  4  He 

F45.  N0+  4  0~*H20  4  He  -  NO  4  H20  4  0  4  He 

F46.  N0+  4  COg^HgO  4  He  ^  NO  4  H20  4  C02  4  0  4  He 

F47.  N0+  4  CO^HgO  4  He  -»  NO  4  Er,0  4  C02  4  02  4  He 
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F17.  0+«H20  +  0“  +  N2  -  02  +  H20  +  0  +  N„ 

F18.  0+*H20  4-  0”  +  N2  -  202  +  H20  4-  N2 
F19.  0^*H20  +  0J  4-  N2  -4  02  +  03  +  H20  +  N2 
F20.  0^*H20  4-  CO"  4-  N2  -*  02  +  0  4-  H20  +  C02  +  N2 

F21.  0j*H20  +  COJ  4  N2  -4  i!02  +  H20  +  C02  +  N2 

F22.  0^*H20  +  02«H20  +  N2  ->  202  +  2H20  +  N0 

F48.  0+«H20  +  0"  +  He  02  +  H20  +  0  +  He 

F49.  o£*H20  +  02  +  He  •*  202  +  H20  +  He 
F50.  OpE20  +  0"  +  He  -*  09  +  0^  +  H20  +  He 
F51.  0  ^H20  +  CO"  +  He  -  02  +  O  +  H20  +  C0<2  4  He 
F52.  o£«H20  +  CO^  +  He  -  202  +  H20  +  C09  +  He 
F53.  0';5«H20  +  0J»H20  +  He  -♦  202  +  2H20  +  He 
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F23.  H304'«H20  4-  O"  4-  N2  -  2H20  +  OH  +  N2 

F24.  H30+«H20  4-  0“  4-  N2  -»  2H20  +  OH  +  0  4-  N2 

F25.  H30'f*H20  4-  O"  -t  N2  -  2H20  4-  OH  4-  02  +  N2 

F26.  H30+*H20  4-  0~  4-  N2  2H20  +  CH  4-  Og  4-  N2 

F27.  H30+«H20  +  C03  +  N2  -  2K20  +  OH  +  C02  +  N2 

F28.  H304'»H20  +  C04  +  N2  -»  2II20  4-  OH  4-  C02  +  O  + 

F29.  H30+«H20  4-  0~*H2Q  +  N?  -  3H20  +  OH  +  0  4-  N'2 
F30.  H30+«H20  +  C0~*H20  +  N2  -  3II20  4-  OH  4-  C02  4-  N2 

F31.  H30+»H20  4-  C0”«*H20  4-  N2  -c  3H20  4*  OH  +  C02  4  O  4  N2 

F54.  H30+*H20  4-  0~  4-  He  -»  2H20  4  OH  4  He 

F55.  H30+»H20  4-  0~  4-  He  -  2^0  4  CH  4  0  4  He 

F56.  H30+»H20  -F  0*3  4-  He  -  2H20  4  OH  4  02  4  He 

F57.  H30+*H20  4-  0~  4-  He  -  2K20  4-  OH  4-  03  4-  He 

F58.  H30+*H20  4-  CO3  4-  He  -»  2H20  +  OH  +  C02  4-  He 

F59.  H30+*H20  4-  C04  4-  He  -  2H20  4-  OH  4-  C02  4  O  4  He 

F60.  H304  •H20  4-  0o»H20  4-  He  -*  3H20  4-  OH  4-  0  4-  He 

F61.  HgO'^HgO  4-  C03*II20  4-  He  -»  3Ho0  4-  OH  4-  C02  +  He 

F62.  H30 f *H20  +  C0'4*H20  4-  He  -  3H<20  4-  OH  4-  C02  4  O  4  He 
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Appendix  B 


PLASMA  LIFETIME  DATA 


The  data  in  Tables  B-l  through  B-4  correspond  to  the  plasma  lifetime 
of  helium  with  air  as  an  impurity  and  of  air  from  sea  level  to  300,000  ft. 
These  data  are  plotted  in  Figures  2  -  5  of  Reference  14  (copy  attached).  The 
numerical  values  in  these  tables  provide  backup  and  constitute  a  database  to 
validate  the  installation  of  AIR.EXE. 

Table  B-l 

PLASMA  LIFETIME  IN  HELIUM  WITH  100  ppm  AIR 


Pressure 

(torr) 

Lifetime  versus  Electron  Number  Density 

lOis/cc 

10“ 

1013 

1012 

ion 

1000 

617  ps 

39.8  n3 

736  ns 

7.4  ps 

25.1  ps 

760 

614  ps 

39.8  ns 

926  ns 

9.3  ps 

29.3  ps 

500 

613  ps 

34.1  ns 

1.2  ps 

12.6  ps 

36.9  ps 

200 

611  ps 

29.3  ns 

2.5  ps 

23.3  ps 

68.1  ps 

100 

611  ps 

27.1  ns 

2.5  ps 

36.9  ps 

117  ps 

50 

- 

27.1  ns 

2.9  ps 

58.4  ps 

200  ps 

20 

- 

27.1  ns 

3.2  ps 

108  ps 

369  ps 

10 

— 

27.1  ns 

3.2  ps 

159  ps 

631  ps 

Pressure 

Lifetime  versus  Electron  Number  Density 

\  lUU  ) 

1010/cc 

109 

108 

107 

106 

117  ps 

342  ps 

584  ps 

631  ps 

631  ps 

1  760 

136  ps 

464  ps 

926  ps 

1.1  ms 

1.1  ms 

159  ps 

736  ps 

1.8  ps 

2.5  ms 

2.5  ms 

251  ps 

1.3  ms 

5.4  pe 

12.6  ms 

15.8  ms 

100 

398  ps 

1.7  ms 

10.0  ps 

31.6  ms 

58.4  ms 

50 

631  ps 

2.5  ms 

13.6  ps 

68.1  ms 

171  ms 

20 

1.1  ms 

5.0  ms 

20.0  ps 

126  ms 

541  ms 

10 

1.6  ms 

7.9  ms 

31.6  ps 

159  ms 

926  ms 

Note:  A  dash  denotes  heating  in  excess  of  50*  C. 
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Table  B-2 

PLASMA  LIFETIME  IN  HELIUM  WITH  10  ppm  AIR 


Pressure 

(torr) 

Lifetime  versus  Electron  Number  Density 

lOis/cc 

1014 

1013 

1012 

ion 

mm, 

616  ps 

42.3  ns 

1.8  ps 

8.6  ps 

73.6  ps 

614  ps 

39.3  ns 

2.1  ps 

10.8  ps 

92.6  ps 

■llfl 

612  ps 

33.9  ns 

2.6  ps 

14.6  ps 

126  ps 

611  ps 

29.0  ns 

3.4  ps 

26.1  ps 

230  ps 

100 

611  ps 

28.7  ns 

3.9  ps 

39.2  ps 

313  ps 

50 

- 

28.6  ns 

3.6  ps 

61.5  ps 

271  ps 

20 

— 

28.7  ns 

3.4  ps 

113  ps 

365  ps 

10 

— 

28.7  ns 

3.3  ps 

178  ps 

627  ps 

Pressure 

Lifetime  versus  Electron  Number  Density 

^  turr  j 

1010/cc 

109 

108 

lO? 

106 

WEm 

mm 

1.2  ms 

3.4  ms 

5.8  ms 

6.3  ms 

HE 

293  ps 

1.4  ms 

4.6  ms 

9.3  ms 

10.8  ms 

369  ps 

1.6  ms 

7.4  ms 

18.5  ms 

25.1  ms 

■1H 

681  ps 

2.3  ms 

12.6  ms 

54.1  ms 

126  ms 

100 

1.2  ms 

3.7  ms 

15.8  ms 

92.6  ms 

316  ms 

50 

1.9  ms 

6.4  ms 

23.3  ms 

126  ms 

631  ms 

20 

3.4  ms 

10.8  ms 

50.1  ms 

171  ms 

1.1  s 

10 

5.4  ms 

15.4  ms 

80.8  ms 

316  ms 

1.4  s 

Note:  A  dash  denotes  heating  in  excess  of  50*  C. 
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Table  B-3 

PLASMA  LIFETIME  IN  HELIUM  WITH  1  ppm  AIR 


Pressure 

(torr) 

■ 

Lifetime  versus  Electron  Number  Density 

1915/cc 

1014 

10*3 

10*2 

10*i 

■EM 

wm 

ps 

42.3 

ns 

2.9 

/iS 

27.0 

(JS 

89.3  (is 

WmM 

614 

ps 

39.3 

ns 

3.1 

/IS 

33.4 

(JS 

112  (is 

mmm 

ps 

33.9 

ns 

3.6 

fJS 

45.1 

(JS 

152  (js 

611 

ps 

29.0 

ns 

4.2 

/JS 

82.5 

(JS 

267  (js 

100 

ps 

28.7 

ns 

4.2 

/ss 

116 

(JS 

414  (js 

50 

— 

28.6 

ns 

3.7 

(JS 

152 

(JS 

655  (js 

20 

- 

28.7 

ns 

3.4 

(JS 

184 

(JS 

1.2  ms 

10 

— 

28.6 

ns 

3.3 

(JS 

226 

(IS 

1.8  ms 

Pressure 

Lifetime  versus  Electron 

Number  Density 

^lun ) 

10lo/cc 

10» 

10» 

107 

106 

mmm 

760 

fJS 

2.5 

ms 

11.3 

ms 

34.2 

ms 

58.4  ms 

■fi m* 

917 

tiS 

2.9 

ms 

13.2 

ms 

47.3 

ms 

92.6  ms 

1.3 

ms 

3.6 

ms 

15.8 

ms 

71.0 

ms 

175  ms 

■1^ 

2.3 

ms 

6.9 

ms 

23.7 

ms 

124 

ms 

540  ms 

100 

2.9 

ms 

11.6 

ms 

37.5 

ms 

158 

ms 

923  ms 

50 

2.5 

ms 

18.4 

ms 

63.8 

ms 

229 

ms 

1.3  s 

20 

3.1 

ms 

32.9 

ms 

107 

ms 

498 

ms 

1.8  s 

10 

4.6 

ms 

46.3 

ms 

152 

ms 

809 

ms 

3.1  s 

Note:  A  dash  denotes  heating  in  excess  of  50*  C. 
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Table  B-4 

PLASMA  LIFETIME  IN  AIR 


Altitude 


Lifetime  versus  Electron  Number  Density 


[Kll) 

— 

1015/cc 

1014 

1013 

1012 

101* 

0 

613 

ps 

4.3 

ns 

8.6 

ns 

9.3 

9.3 

ns 

10 

651 

ps 

7.4 

ns 

17.1 

ns 

20.0 

20.0 

ns 

20 

662 

ps 

11.7 

ns 

31.6 

ns 

i m 

30 

668 

ps 

17.1 

ns 

63.1 

ns 

H 

40 

679 

ps 

21.5 

ns 

117 

ns 

171 

ns 

200 

ns 

50 

693 

ps 

23.3 

ns 

215 

ns 

398 

ns 

501 

ns 

60 

702 

ps 

23.3 

ns 

369 

ns 

858 

ns 

1.3 

ps 

70 

- 

23.3 

ns 

501 

ns 

1.7 

/iS 

2.9 

ps 

80 

- 

25.1 

ns 

631 

ns 

3.2 

ps 

6.8 

ps 

90 

— 

27.1 

ns 

681 

ns 

5.4 

ps 

13.6 

ps 

100 

— 

29.3 

ns 

631 

ns 

7.4 

ps 

25.1 

ps 

— 

— 

29.3 

ns 

631 

ns 

8.6 

ps 

43.0 

ps 

Wmm-- 

- 

— 

631 

ns 

8.6 

ps 

63.1 

ps 

- 

— 

681 

ns 

8.6 

ps 

73.6 

ps 

140 

- 

- 

736 

ns 

7.9 

ps 

85.8 

ps 

150 

- 

- 

736 

ns 

7.9 

IIS 

85.8 

ps 

160 

- 

— 

736 

ns 

7.9 

US 

85.8 

ps 

170 

- 

— 

794 

ns 

7.9 

1* 

85.8 

ps 

180 

- 

•- 

— 

7.9 

ps 

79.4 

ps 

190 

- 

— 

- 

8.6 

m 

79.4 

ps 

200 

— 

— 

8.6 

ps 

79.4 

ps 

a— 

_ 

— 

— 

9.3 

ps 

79.4 

ps 

- 

- 

-- 

9.3 

ps 

S5.8 

ps 

- 

- 

- 

9.3 

ps 

85.8 

ps 

- 

- 

— 

- 

85.8 

ps 

250 

- 

— 

— 

— 

85.8 

ps 

IlEi&iH 

- 

— 

— 

- 

92.6 

ps 

Kml  ■ 

- 

— 

— 

— 

92.6 

ps 

— 

-- 

— 

- 

— 

- 

— 

— 

— 

■  K 

— 

— 

— 

— 

— 

Note:  A  dash  denotes  heating  in  excess  of  50*  C. 


Table  B-4  (Concluded) 


Altitude 

(kft) 


Lifetime  versus 

Electron  Number  Density 

1010/cc 

109 

108 

107 

106 

9.3 

ns 

20.0 

ns 

43.0 

ns 

85.8 

ns 

200 

ns 

501 

ns 

1.4 

pis 

3.4 

pa 

8.6 

pa 

21.5 

pa 

46.4 

pa 

92.6 

pa 

171 

271 

pa 

369 

pa 

464 

pa 

584 

pa 

736 

pa 

926 

pa 

926 

pa 

858 

pa 

794 

pa 

794 

pa 

794 

pa 

794 

pa 

794 

1* s 

858 

pa 

858 

pa  j 

858 

pa 

858 

pa 

858 

ifj 

0.3  ns 
20.0  ns 
43.0  ns 
85.8  ns 
200  ns 
501  ds 

1.4  pa 

3.4  pa 

9.3  pa 

23.3  pa 

58.4  us 


858  pa 
1.3  ms 


9.3 

ms 

9.3 

m3 

8.6 

ms 

8.6 

ms 

7.9 

ms 

7.9 

ms 

7.9 

ms 

7.9 

ms 

7.9 

ms 

7.9 

ms 

9.3  n3 
20.0  ns 
43.0  ns 
85.8  ns 

200  ns 
501  ns 

1.4  pa 

3.4  pa 

9.3  pa 

23.3  pa 
63.1  pa 

159  pa 
398  ps 
858  pa 
1.7  ins 


100  ms 
100  ms 


200 

ns 

501 

ns 

1.4 

pa 

3.4 

pa 

9.3 

pa 

23.3 

pa 

63.1 

ps 

159 

pa 

398 

pa 

1.0 

ms 

2.3 

ms 

5.0 

ms 

9.3 

ms 

17.1 

ms 

31.6 

ms 

58.4 

ms 

oo 

o 

ms 

200 

ms 

342 

ms 

584 

ins 

794 

ms 

926 

ms 

1.0 

s 

1.0 

s 

926 

ms 

858 

ms 

794 

ms 

200  ns 
501  ns 

1.4 

3.4 
9.3 

23.3 
63.1 


159 
398 
1.0 

2.3 

5.4 
11.7 

23.3 

46.4 
100 
184 


341  ms 
681  ms 


Appendix  C 


AIR  CHEMISTRY  CODE  ORGANIZATION 

The  air  chemistry  code  has  been  developed  to  run  on  an  IBM  AT 
computer  equipped  with  a  math  coprocessor  chip.  The  source  code  was 
developed  using  IBM  FORTRAN,  Ryan-McFarland  Corp.  version  1.22.  The 
software  provided  as  an  attachment  to  this  report  includes  source  code,  EXE 
files,  and  auxilary  programs.  It  is  on  a  high-density,  5f-in.,  96-TPI,  1.2-MB 
diskette  and  requires  a  1.2-MB  high-density  disk  drive. 

The  source  code  is  provided  to  other  researchers  via  AFOSR/NP  on  an 
as-is  basis;  that  is,  no  guarantee,  no  warranty,  no  software  support,  and  no 
hardware/installation  support . 

If  other  equipment  is  used  or  if  another  FORTRAN  compiler  is  used 
some  source  code  changes  will  be  necessary.  For  example,  Ryan-McFarland 
Corp.  version  1.00  (1984)  had  several  integer  and  real*4  to  real*8  conversion 
bugs  that  severely  affected  results  from  the  air  chemistry  code  for  altitudes 
above  200,000  ft.  The  difficulty  was  not  in  overflow  (+++++)  or  underflow 

( - )  5ut  jn  not-a-number  (NAN)  generation  (?????).  Version  1.22  and 

better  source  code  programing  have  eliminated  this  problem.  The  data  in 
Appendix  B  for  helium  and  air  should  be  duplicated  before  any  modifications 
to  the  source  code  are  attempted. 

Variable  names  were  chosen  to  allow  easy  interpretation  of  the  source 
code.  Comments  are  included  to  help  simplify  code  modification.  Table  C-I 
is  a  succinct  guide  to  how  tne  program  was  organized  and  manipulated  into 
64-kB  blocks,  t  hree  common  blocks  were  coded  (as  include  files)  to  simplify 
the  exchange  of  data  between  subroutines.  Coding  these  common  blocks  as  a 
matrix  and  passing  them  as  an  explicit  variable  was  initially  attempted.  This 
approach  eventually  produces  an  error  related  to  the  number  of  variable  names 
in  a  subroutine  call.  The  large  common  blocks  eliminated  that  difficulty. 

The  function  and  its  Jacobian  had  another  subroutine  limitation:  the  number 
of  instructions.  This  was  overcome  by  breaking  large  subroutines  into  smaller 
ones,  and  then  calling  them  sequentially. 
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Table  C-l 

PROGRAM  ORGANIZATION 


DOS  File  Name 

Subroutine 

Common  Blocks 

AIRA.FOR 

-  include  - 

/AIRSMOL/ 

AIRB.FOR 

-  include  - 

/RATES/ 

AIRC.FOR 

-  include  - 

/EXTENT/ 

AIR.  FOR 

Main  program 

-  none  - 

TABLESSA 

/blank  common/, /AIRSMOL/ 

ENG 

-  none  - 

AIRI.FOR 

PHOTO 

/AIRSMOL/ 

AIRS^OMP 

/AIRSMOL/ 

XRA 

/AIRSMOL/ 

MOMENTUM 

/AIRSMOL/ 

FSRATE 

-  none  - 

SETSRATE 

-  none  -  (ENTRY  in  FSRATE) 

AIR2.FOR 

ASRATES 

/AIRSMOL/, /RATES/ 

AIR3.FOR 

AIRSCHEM 

/blank  common/, /AIRSMOL/ 

AIR4.FOR 

FSODE 

FSYDOT 

/AIRSMOL/, /RATES/, /EXTENT/ 
/AIRSMOL/, /EXTENT/ 

AIR5.FOR 

JSODE 

/AIRSMOL/, /RATES/ 

J$A 

/AIRSMOL/, /RATES/ 

JSB 

/AIRSMOL/, /RATES/ 

JSC 

/AIRSMOL/, /RATES/ 

AIR6.FOR 

JSNEG 

/AIRSMOL/, /RATES/ 

JSPOS 

/AIRSMOL/,/ RATES/ 

CBALANCE 

-  none- 

ODE1.FOR 

-  include  - 

/AIRSMOL/ 

ODE2.FOR 

-  include  - 

/RATES/ 

ODE3.FOR 

-  include  - 

/EXTENT/ 
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Source  code  for  the  air  chemistry  program  is  divided  into  the  13  files 
described  in  Table  C-l.  The  description  of  each  program  element  in  Table 
C~2  will  help  in  any  code  modification.  Modification  of  this  code  is  complex, 
even  after  the  operation  of  LSODE  is  understood.  Addition  or  reduction  of  a 
species  involves  a  change  in  the  number  of  unknowns,  which  affects  the 
dimensioning  of  y,  ydot,  and  pd;  subroutines  CBALANCE  and  TABLESSA; 
the  output  data  matrix;  and  the  AIRSMOL,  RATES,  and  EXTENT  common 
blocks.  Adding  to  or  reducing  the  number  of  reactions  is  much  simpler, 
because  only  y,  ydot,  pd,  ASRATES,  and  the  EXTENT  common  block  are 
affected . 

Programming  can  be  profoundly  simplified  by  using  the  numerical 
Jacobian  internal  to  LSODE;  however,  the  execution  time  of  a  simulation 
increases  significantly.  This  possibility  is  only  suggested  for  operation  on  a 
machine  that  is  at  least  an  order  of  magnitude  faster  than  an  early  model 
IBM  AT.  The  explicit  Jacobian  coding  was  necessary  to  reduce  the  execution 
time  for  a  simulation  by  a  factor  of  50.  The  data  in  Appendix  B  took 
approximately  5  min  per  data  point. 

Auxiliary  software  on  the  diskette  is  noted  in  Table  C-3.  The 
differential  equation  solver  used  in  this  program  is  LSODE,  developed  by 
Hindmarsh  [34].  A  copy  of  the  support  documentation  for  the  Hindmarsh  code 
is  provided  in  ODEO.TXT.  The  program  JACOBIAN  checks  the  coding  of  the 
explicit  Jacobian.  If  any  changes  are  made  to  the  number  of  reactions  or  the 
number  of  species,  JACOBIAN  should  be  run  to  verify  complete  coding. 

JACOBIAN  was  developed  because  some  simulation  data  did  not  seem  to 
make  sense.  A  number  of  species  were  continually  assuming  negative  values 
when  there  was  no  physical  reason  for  it.  The  difficulty  was  traced  back  to 
eight  coding  errors  in  1500  lines  of  the  Jacobian  and  a  common  block  that 
was  coded  differently  in  one  subroutine.  The  include  statements  were  adopted 
to  eliminate  the  latter  difficulty.  LINKSJAC.BAT  is  an  example  of  how  to 
link  all  the  required  subroutines  of  JACOBIAN. 

Because  the  program  contains  many  lines  of  code,  each  of  13  files  can  be 
compiled  separately  to  save  time.  LINKSAIR.BAT  is  an  example  of  how  to 
link  these  separate  segments  into  an  EXE  file. 
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Table  C-2 

PROGRAM  ELEMENTS 


Program  Element  Description 

Main  program  Sequentially  calls  subroutines  to  define  gas 

composition,  initial  conditions,  and 
differential  equation  solver  LSODE 

Output  FORT70  Written  to  default  disk  drive,  contains 

tabular  output  of  simulation  identical  to  the 
output  directed  to  the  screen 

Output  FORT71  Written  to  default  disk  drive,  contains 

numerical  simulation  data  saved  in  the  format 
of  I,AIR$DATA(60)  described  in  AIRSCHEM 

/AIRSMOL/  Common  block:  concentration  of  gas  species 

/RATES/  Common  block:  reaction  rates  keyed  to  numbering 

system  in  Appendix  A 

/EXTENT/  Common  block:  extent  of  reaction  with  numbering 

keyed  to  reactions  in  Appendix  A 

TABLESSA  Subroutine:  Output  data  matrix  to  screen  and  to 

printer  via  the  output  file  FORT70 

ENG  Converts  Real *4  numbers  to  engineering  notation 

PHOTO  Photoionization  data  for  atmospheric  ions  and 

a  xenon  flashlamp  (incomplete) 

AIRSCOMP  Initializes  neutral  gas  densities  for  earth's 

atmosphere,  synthetic  air,  and  helium 

XRAY  Initializes  ionization  by  product  densities 

data  extracted  from  various  references 

MOMENTUM  Momentum-transfer  collision  rate  based  on  pas 

composition  and  temperature,  Itakawa  [18] 

FSRATE  Sets  parameters  that  define  an  ionization  pulse: 

ionization  rate,  duration,  waveform,  etc. 

SETSRATE  ENTRY  in  FSRATE  for  pulsed  ionization 

A$ RATES  Defines  reaction  rates 

AIRSCHEM  Defines  Y  and  AIRSDA  TA  matrices 

Selects  simulation  parameters 
CALL  to  LSODE  solver 
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Table  C-2  (Concluded) 


Program  Element 


ODE3 


Description 


Subroutine:  contains  the  extent  of  reaction 
and  a  CALL  to  FSYDQT  that  defines  the 
derivatives  of  the  function  that  LSODE  solves 
FSODE  is  defined  as  an  EXTERNAL 
The  name  FSODE  is  passed  in  a  CALL  to  LSODE 
in  subroutine  AIRSCHEM  (AIR3.FOR) 

FSODE  is  passed  to  LSODE  in  subroutine 
AIRSCHEM  in  AIR3.FOR 
Y  is  the  species  concentration  matrix 
YDOT  is  the  differential  equation  matrix 


YDOT  derivatives  related  to  extent  of  reaction 


Jacobian  call  by  LSODE:  method  flag  21 

User-supplied  full  Jacobian 

JSODE  is  defined  as  an  EXTERNAL, 

The  name  JSODE  is  passed  in  a  CALL  to  LSODE 
in  coroutine  AIRSCHEM  (AIR3.FOR) 

JSODE  consists  of  seven  segments 
JSODE  contains  part  of  the  neutral  species 


Part  of  JSODE,  neutral  species 


Part  of  JSODE,  neutral  species 


Part  of  JSODE,  neutral  species 


Part  of  JSODE,  negative  species 


Part  of  JSODE,  positive  species 


Part  of  JSODE,  modifies  values  of  Y  mat^x 
to  maintain  charge  equilibrium 


Lawrence  Livermore:  ordinary  differential 
equation  solver  (Hindmarsh);  documentation, 
control  codes,  and  error  codes  on  file 
ODEO.TXT,  first  of  three  segments 


LSODE:  second  of  three  segments 


LSODE:  third  of  three  segments 
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Table  C-3 

AUXILIARY  SOFTWAFE 


DOS  File  Name 


ODEO.TXT 


LINKSAIR.BAT 


AIR.EXE 


JACOBIAN.FOR 


LINKSJAC.BAT 


JACOBIAN.EXE 


CCDR.FOR 


Description 


Documentation  for  LSODE:  control  codes,  error 
codes,  example  program 


BAT  file  to  link  individual  program  segments 
into  AIR  EXE 


EXE  file  for  air-chemistm  code 


An  explicit  Jacobian  checker  to  check  accuracy 


BAT  file  to  link  individual  program  segments 
into  JACOBIAN.EXE 


EXE  file  for  JACOBIAN 


CCDR.EXE 


EPSTEIN.FOR 


Cold-Collisional-Dispersion  relation  versus 
frequency  versus  electron  density  versus 
collis’on  frequency 


EXE  file  for  CCDR 


EPSTEIN.EXE 


Scattering  and  absorption  by  an  Epstein  profile 
Generates  ASCII  file  FORT70:  reflection 
coefficient  or  total  return  [14]  versus 
frequency 


EXE  file  for  EPSTEIN 


CCDR  is  provided  as  a  means  of  evaluating  tire  cold-collisional- 
dispersion  relation  for  any  frequency,  collision  rate,  and  electron  number 
density.  It  is  easy  to  run  and  provides  tabular  listings  in  file  FORT70  of  the 
complex  wave  number  as  a  function  cf  the  above  quantities.  This  program 
can  be  run  for  a  single  frequency  and  collision  rate  with  the  number  density 
set  to  vary.  The  resulting  variation  in  wave  number  as  a  function  of  number 
density  can  be  interpreted  as  the  effect  on  wave  number  of  a  plane  wave 
incident  on  an  Epstein  profile.  If  Re( k)  remains  close  to  its  free-space  value, 
the  backscatter  reflection  coefficient  is  small.  If  /m(k)  increases,  the  plasma 
will  exhibit  some  attenuation. 
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EPSTEIN  is  a  program  that  computes  the  reflection  or  round-trip 
absorption  of  a  plane  wave  on  an  Epstein  profile  as  a  function  of  frequency. 

It  is  an  implementation  of  the  reflection  coefficient  and  total  reduction 
described  in  Section  IV  of  Reference  14  (copy  attached).  This  program 
prompts  for  reflection  or  total  reduction,  collision  rate,  frequency  range, 
electron  number  density,  and  the  Epstein  profile  gradient  scale  factor.  Up  to 
four  sets  of  electron  densities  and  gradient  scale  factors  can  be  described  each 
time  this  program  is  run.  The  output  of  EPSTEIN  is  ASCII  file  FORT70, 
which  is  in  a  format  suitable  for  most  plotting  packages. 
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Abstract  —  Tenuous  plasmas  in  the  earth's  atmosphere  from  sea  level  to 
100  km  can  be  modeled  as  a  cold  collisionai  plasma.  If  an  ionization  process 
results  in  an  electron  density  profile,  such  as  an  Epstein  profile,  that  decreases 
from  a  maximum  as  a  function  of  distance,  then  the  plasma  can  either  reflect 
or  absorb  electromagnetic  waves,  depending  on  the  plasma  characteristics.  A 
high  reflection  coefficient  requires  a  grazing  angle  of  incidence,  a  low 
momentum-transfer  collision  rate,  a  high  plasma  density,  and  a  short  free- 
space  to  plasma  transition  relative  to  a  wavelength.  High  absorption  from 
VHF  through  \-band  requires  a  high  collision  rate,  a  low  plasma  density,  and 
a  plasma  transition  of  approximately  one  wavelength.  Both  extremes  are 
discussed.  Typical  collision  rates  and  plasma  lifetimes  at  atmospheric  pressure 
are  quantified.  Power  required  to  sustain  a  plasma  in  air  is  high  because  of  a 
short  plasma  lifetime,  but  a  noble  gas  contained  within  a  membrane  has  a 
longer  plasma  lifetime  and  a  lower  power  requirement.  Two  means  of 
ionization  are  discussed:  high-energy  electron-beam  impact  and  UV 
photoionization  of  an  organic  vapor. 

Supported  by  Air  Force  Office  of  Scientific  Research,  AFOSR/NP, 
contract  F49620-85-K-0013. 


1 


I.  INTRODUCTION 


The  electromagnetic  properties  of  a  plasma  at  atmospheric  pressure  is  an 
interdisciplinary  topic  that  combines  plasma  physics,  radio  wave  propagation, 
and  air  chemistry.  In  this  paper,  the  atmosphere  extends  from  sea  level  to 
100  km  (300,000  ft),  and  the  plasma  body  exists  in  free  space.  Prior  research 
on  naturally  occurring  and  man-made  atmospheric  plasma  bodies  quantifies 
radio  wave  effects.  Gerson  [1]  describes  EMayer  absorption,  Gunar  and 
Mennella  [2]  note  fluctuations  in  radar  backscatter  as  well  as  the 
communication  blackout  of  reentry  vehicles,  and  Glasstone  and  Dolan  [3] 
describe  the  effects  of  a  nuclear  weapon  generated  plasma  on  radio  and  radar 
signals. 

Figure  1  depicts  two  uses  of  man-made  plasmas.  Borisov  and  Gurevich 
[4]  and  Gurevich  [5]  suggested  that  crossed  high-power  electromagnetic  pulses 
could  break  down  air  below  the  D-layei  and  so  form  a  plasma  layer.  This 
layer  could  be  used  to  reflect  signals  to  a  great  range,  and  the  highest 
frequency  for  reflection  would  depend  on  the  electron  number  density  in  the 
slab.  A  plasma  could  also  be  used  as  an  absorber,  as  illustrated  in  Fig.  1(b). 
The  absorption,  power  required,  and  bandwidth  are  described  in  this  paper. 

These  man-made  plasmas  each  use  an  ionization  source  to  generate  a 
tenuous  plasma,  resulting  in  an  electron  number  density,  n  ,  that  is  high  near 
the  source  and  diminishes  with  distance  from  the  source.  After  the  source 
shuts  off,  nQ  decreases  as  a  function  of  time  as  electrons  recombine  with 
positive  ions  or  attach  to  form  negative  ions.  Section  II  describes  the 
electromagnetic  properties  that  are  essential  to  an  understanding  of  these 
man-made  plasmas.  The  theory  model  for  scattering  parameters  is  derived 
from  the  work  of  Epstein  [6]  as  presented  by  Budden  (7)  in  the  context  of 
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scattering  from  plasma  gradients  in  the  ionosphere.  Electromagnetic  properties, 
such  as  wave  numbei ,  power  reflection  coefficient,  and  absorption  coefficient, 
are  a  function  of  the  host  gas  and  plasma. 

Section  III  presents  gas  and  plasma  characteristics,  such  as  momentum- 
transfer  collision  rate,  plasma  lifetime,  recombination  kinetics,  and  the  effect  of 
noble  gases.  For  example,  a  noble  gas,  used  instead  of  air  as  a  host  gas, 
extends  the  plasma  lifetime  by  reducing  negative-ion  formation.  Several  tables 
and  figures  quantify  a  wide  range  of  values  used  to  estimate  the  performance 
and  power  required  to  sustain  a  plasma. 

Applications  for  a  plasma  with  a  gradient,  discussed  in 
include  (a)  a  high-altitude  plasma  that  can  reflect  or  absorb  from  HF  to 
VHF,  and  (b)  a  broadband  atmospheric  pressure  absorber.  A  novel 
characteristic  of  plasma  used  as  an  absorber  is  its  bandwidth  that  extends 
from  VHF  through  X-band  at  atmospheric  pressure. 

Section  V  describes  the  generation  and  use  of  plasma,  including  electron 
impact  ionization  with  a  high-energy  electron-beam  source  and  UV 
photoionization  of  an  alkali  vapor  or  an  organic  vapor,  such  as 
tetrakis(dimethylamino)ethylene,  TMAE,  [8}— {10].  The  power  required  to 
sustain  a  plasma  is  quantified,  and  properties  such  as  maximum  absorption 
and  bandwidth  are  discussed.  Trade-offs  among  maximum  absorption, 
absorption  bandwidth,  duty  ratio,  and  applied  power  permit  optimization  of 
absorption  primarily  at  VHF. 


II.  ELECTROMAGNETIC  CHARACTERISTICS 


A.  Cold  Collisional  Plasma 

In  an  atmospheric  pressure  plasma,  electrons  undergo  numerous  <^cllisions 
with  atoms  and  molecules  in  the  neutral  background  gas.  Collisions  damp 
electron  motion  and  so  convert  electromagnetic  energy  coupled  into  electron 
motion  directly  to  heat.  Electrons  produced  by  an  ionization  source  cool 
rapidly  in  collisions  with  the  background  gas.  The  host  gas  (air  or  a  noble 
gas  confined  by  a  membrane)  is  at.  the  ambient  atmospheric  temperature  and 
under  many  circumstances  has  sufficient  heat  capacity  to  cool  electrons  with 
little  change  in  temperature. 

If  the  ionization  process  maintains  an  electron  temperature  above  1  eV, 
then  a  multi-specie  full-kinetic  model  is  necessary  because  of  electronic 
excitation,  dissociation,  ionization,  an  so  on.  Whenever  host  gas  heating  due 
to  generation  of  a  tenuous  plasma  body  is  negligible,  then  the  electron 
temperature  is  near  the  ambient  temper&ti^e.  After  hot  electrons  cool  to  a 
fraction  of  an  eV,  plasma  characteristics  can  be  approximated  via  a  cold 
collisional  plasma  model  [llj,  [12]. 

The  dispersion  relation  for  a  cold  collisional  plasma  with  a  constant 
momentum-transfer  collision  rate,  v,  is 


k(w) 


u(u  -  \u)j 


l/? 


(1) 


where  k  is  a  complex  wave  number,  v  is  angular  frequency,  i  =  JT,  and  c  is 

2  1 12 

the  speed  of  light.  The  plasma  angular  frequency  is  =  (n  e  /m  e  )  '  , 

V  V  U 

where  e  and  m  are  the  electron  charge  and  mass,  n  is  the  electron  number 
density,  and  eQ  is  the  free-space  permittivity. 


The  wave  number  is  k  =  kf  -f  ikj,  where  kr  —  /2e(k)  and  kj  =  /m(k). 

Wave  motion  is  proportional  to  exp[-ri(wt-krz)j  exp(-fkjz),  where  z  is 

propagation  distance  and  kj  is  the  attenuation  constant.  A  lossless  plasma 

(1/  =  0)  exhibits  a  wave  propagation  cutoff  for  u  <  ay  A  collisional  plasma 

with  n  =  2  *  10^  cm  ^  and  v  =  300  *  10^  s  *  has  w  =  25  *  10^  s  *,  a 
e  p 

conductivity  similar  to  sea  water,  and  a  wave  number  that  remains  finite  for 

u>  <  u)  .  The  value  of  k  at  1  GHz  is  k  =  20.96  and  k.  =  -3.475  m~V 
p  r  i 
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For  u  <  v  and  u  <  oy/i/,  kr  and  kj  simplify  to  kf  ~  u/c  and 

o 

kj  2  -a>p/(2 vc).  Absorption  to  first  order  is  broadband  in  nature  and 
independent  of  frequency  below  u  jv.  But,  the  transition  from  free  space  to 
plasma  produces  a  reflection  that  is  frequency  dependent.  The  magnitude  of 
the  reflection  depends  on  the  spatial  distribution  of  plasma  in  the  transition 
region. 


B.  Epstein  Profile  Model 

An  Epstein  profile  provides  an  analytical  model  for  estimating  scattering 
coefficients  from  plasmas  tbai  have  a  smooth  electron  density  gradient  [6],  [7j. 
An  Epstein  electron  n  imber-density  profile  is 

o's)  - - ^ -  (2) 

1  +  exp  (— z/ cr) 

where  n0  is  the  electron  density  at  z  =  +  ®,  y  is  a  gradient  scale  factor,  and 
v  is  independent  of  z.  A  plane  wave  propagating  from  z  =  -®  to  z  =  +® 
encounters  an  increasing  electron  number  density  that  results  in  a  reflection 
determined  by  the  gradient  scale  factor,  a,  and  nQ. 


C.  Power  Reflection  and  Transmission  Coefficients 


The  power  reflection  and  transmission  coefficients  for  a  wave  incident  at 
angle  6  on  an  Epstein  profile  are  described  by  Budden  [7,  Chapter  11] 


C  -q 

2 

r[l  4-  ik<r(q  4-  C)] 

C  +  q 

T[1  4-  ik<x(q  —  C)] 

(3) 


T  = 


|Cq| 


|C  4-  q|2 


f*[l  +  ik<r(q  +  C)j 


r[l  4-2ik<rq]r[l  4-  2ik<7C] 


exp[4-2/m(q)z]  (4) 


2  _  r2 
q  =  C  - 


u(u )  -  iy) 


(5) 


where  k  is  the  free-space  wave  number,  C  =  cos (5),  and  q  is  a  solution  of  the 
Booker  quartic  described  in  [7]. 

The  second  term  of  q2  in  (5)  has  a  magnitude  much  less  than  unity  at 
high  frequency,  w  >  wp  and  u  >  i/,  but  is  large  and  primarily  complex  for 
u>  <  a/p  and  w  <  v.  Before  resorting  to  the  machine  calculations  of  (3)  and 
(4)  in  Section  IV,  consider  three  limiting  cases. 

For  o  -*  0,  an  Epstein  profile  approximates  a  slab  discontinuity  at 
z  =  0.  The  arguments  of  the  gamma  functions  in  (3)  and  (4)  equal  unity, 
and  the  values  of  R  and  T  reduce  to  the  Fresnel  reflection  formulas  [7]. 

The  second  case  is  for  grazing  angle  illumination,  0  -*  90*  so  that  C  -*  0. 
The  arguments  of  the  gamma  functions  in  (3)  and  (4)  then  approach  unity  for 
any  o  and  v.  The  reflection  coefficient  approaches  unity,  and  the  transmission 
coefficient  approaches  zero.  This  suggests  the  use  of  a  tenuous  plasma  as  a 
reflector  at  a  grazing  angle. 
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The  third  case  is  for  0=0,  the  backscatter  reflection  coefficient.  In  (3) 
the  argument  of  the  gamma  function  in  the  numerator  is  4-  2ik<r  larger  than 
that  in  the  denominator.  Using  Euler's  formula,  the  absolute  value  for  the 
ratio  of  gamma  functions  is 


r<»  +  aka)"  _  lim 

z 

4 

(z+1) 

4 

(z+n) 

r(z)  (  nHfl0 

(z+iy) 

(z+  i  y-t-1) 

(z+  i  y+n) 

where  y  =  2k<r  and  z  =  1  +  ikff(q-l).  For  u  <  wp  and  u  <  v,  the  value  of 
q  reduces  to  q  :  C  -  w*/(2u>vC)  or  q  -  1  =  -iu£/(2av)  for  0=0.  For  any 

r  r 

value  of  a  >  0,  z  is  real  and  y  is  imaginary.  Therefore  each  term  in  the 
infinite  product  maintains  |(z+n)/(z-Hy+n)|4  <  1  and  the  value  of  (6)  is  less 
than  one. 

Because  R  is  the  product  of  the  Fresnel  reflection  coefficient  ana  (6), 

reflection  from  an  Epstein  profile  varies  with  frequency  and  slab  thickness, 
o 

For  u  <  u;p/(4i/),  the  wavelength  is  much  longer  than  the  slab  dimension, 
z  >  y,  the  value  of  (6)  is  approximately  one,  and  the  reflection  coefficient 
reduces  to  the  Fresnel  term.  For  oj  =  and  a  =  A/(4x),  y  =1,  z  =  2, 

and  the  first  few  terms  are  propo.  tional  to  (4/5)2.  At  a  higher  frequency  such 
that  u  =  io p/y  and  a  =  A/4,  y  =  x,  z  =  i  +  x/4,  and  the  first  few  terms 
are  proportional  to  "(1/4)  .  Hence,  the  reflection  coefficient  for  backscatter 
can  become  small.  A  small  reflection  coefficient  results  from  incoherent 
backscatter  from  the  plasma  gradient  rather  than  a  coherent  reflection  at  a 
slab  discontinuity. 
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D.  Attenuation  Coefficient 


In  addition  to  scattering  fro  the  plasma,  a  signal  that  propagates  into 
an  Epstein  gradient  and  then  reflects  from  a  g<od  conductor  will  have  its  field 
quantities  reduced  by  exp(-2|kjdz),  where  the  integral  is  from  z  —  -®  to  the 
surface  of  the  conductor  at  +zQ.  The  factor  of  2  in  the  exponential  denotes  a 
round  trip.  For  a  tenuous  plasma,  k.  is  proportional  to  n  integrated  over  all 

I  C 

space.  The  integral  of  an  Epstein  profile  from  — zQ  to  +zQ  is  n0zQ.  The 

2 

round-trip  attenuation  in  dB  for  a  plasma  with  u  <  u  and  u>  <  is 


A(dB)  =  4.343 


'  e2b  ' 

r  a 

je 

t  m  c 

V 

o  e 

j 

where  h  (=2zQ)  is  the  total  thickness  of  the  plasma  and  ng  is  in  m 


-3 


(7) 


The 


5%  and  95%  values  of  n(z)  for  an  Epstein  profile  correspond  to  tzq  -  2.944a, 
respectively.  For  an  Epstein  profile  spanning  the  5%  to  95%  range,  h  is 
5.888a,  and  attenuation  is  A(dB)  =  2.71  x  10-4  n  ajv. 


E.  Numerical  Evaluation 


To  evaluate  the  reflection  and  transmission  coefficients  in  (3)  and  (4), 
the  gau'jma  function  T(z)  for  complex  arguments  in  the  right-half  complex 
plane  must  be  evaluated.  Lanczos  [13]  developed  an  asymptotic  expansion 
correct  to  better  tha^  0.1%  for  the  right-half  complex  plane.  His  formula  is 


f(z+!)  =  (z  +  2)^+1/2^  e"(z+2)  (2*-)^ 


0.999779  + 


1.084635 


z  +  1 


(8) 


Because  the  magnitude  of  z  can  become  large  in  (3)  and  (4),  computing 
ln[r(z+l)]  prevents  an  overflow  condition  and  retains  accuracy. 
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III.  PLASMA  CHARACTERISTICS 


A.  Momentum- Transfer  Collision  Rate 

The  electron  momentum-transfer  collision  rate,  v,  is  a  function  of 
temperature  and  gas  species  and  sets  the  value  of  kj  that  determines 
electromagnetic  absorption.  Itikawa  [14j  discusses  effective  values  of  v  for 
atmospheric  gases  and  noble  gases  over  the  electron  temperature  range  of 
100"  K  to  5000*  K;  the  data  in  Table  I  spans  300*  K  to  1000*  K.  In  air,  20,000 
collisions  vill  reduce  the  electron  energy  by  a  factor  of  40  [11].  pressure 

of  760  tori  (sea  level)  this  process  requires  50  ns,  based  on  the  1000*  K  rate. 

If  thermalization  is  faster  than  or  comparable  to  a  plasma  lifetime,  then  an 
effective  electron  temperature  during  the  afterglow  is  near  the  ambient 
t°rT»erature  of  the  host  gas. 


TABLE  I 

MOMENTUM -TRANSFER  COLLISION  RATE 


Gas  Species 

i//N 

[lO"8  cm3s  l) 

Tg  =  300*  K 

500*  K 

1000* K 

n2 

0.594 

0.959 

1.77 

02 

0.286 

0.437 

0.804 

co2 

10.05 

9.68 

7.48 

h2o 

77.90 

56.25 

34.6 

He 

0.763 

1.01 

1.50 

Ne 

0.076 

0.118 

0.217 

Ar 

0.220 

0.145 

0.094 

Kr 

1.775 

1.315 

0.713 

Xe 

5.29 

3.81 

1.90 

Dry  Air 

0.529 

0.845 

1.550 

Data:  N  is  gas  density  in  cm'3 
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At  a  pressure  of  760  torr  and  an  electron  temperature  of  5Q0CK,  the 

9—1  9—1 

collision  rates  for  dry  air  and  helium  are  212  *  10  s  and  253  *  10  s  . 
Water  vapor  strongly  influences  the  collision  rate  at  low  altitudes.  For 
example,  50%  relative  humidity  increases  the  sea-level  dry-air  collision  rate 
from  212  *  IQ9  s-1  to  417  *  109  s"1.  The  values  of  v  in  Table  II  include 
the  variation  of  air  composition  with  altitude  from  sea  level  to  300,000  ft 
(91.5  km). 


TABLE  II 

COLLISION  RATE  IN  AIR 


Altitude 

(kf* 

Pressure 

(torr) 

Momentum  Transfer,  Collisions  per  sec 

Te  =  300*  K 

500*  K 

1000* K 

0 

760 

466  *  109 

417  x  109 

425  *  109 

11.3 

500 

103 

121 

180 

32.8 

200 

23.6 

37.6 

72.0 

47.5 

100 

12.1 

19.3 

37.1 

62.0 

50 

5.07 

9.70 

18.6 

81.4 

20 

2.38 

3.80 

7.27 

96.4 

10 

1.18 

1.87 

3.57 

111.6 

5 

0.584 

0.914 

1.73 

133.4 

2 

0.272 

0.382 

0.672 

150.5 

1 

0.216 

0.248 

0.365 

169.1 

500  x  10" 

3  122  «  106 

134  x  106 

?89  x  106 

192.4 

200 

26.7 

37.8 

67.1 

209.8 

;oo 

11.8 

18.3 

34.4 

225.5 

50 

5.92 

9.42 

18.0 

245.1 

20 

2.49 

3.98 

7.64 

253.7 

10 

1.27 

2.03 

3.89 

274.1 

5 

0.662 

1.06 

2.03 

291.1 

2 

0.257 

0.411 

0.789 

300.0 

1.32 

0.161 

0.258 

0.494 
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B.  Plasma  Lifetime 

The  plasma  lifetime,  r.  is  the  time  required  for  a  plasma  of  initial 
density  nQ  to  decrease  in  concentration  by  a  factor  of  1/e.  It  is  useful  in 
estimating  the  power  required  to  sustain  a  plasma.  Estimates  of  r  have  been 
computed  with  an  air  chemistry  code  developed  to  model  the  deionization 
process  [15],  [16].  This  code  takes  into  account  the  major  atmospheric  gases 
plus  B^O  as  a  function  of  altitude.  The  same  code  has  an  option  to  treat 
helium  as  the  major  species  with  air  as  an  impurity. 

The  Earth’s  Atmosphere:  The  plasma  lifetime  in  ail  exhibits  a 

dependence  on  electron  number  density  and  on  pressure  (altitude)  as  shown  in 

12  —'1 

Fig.  1.  For  ng  <10  cm  ,  three-body  attachment  to  oxygen  dominates  the 
deionization  process.  The  slope  of  the  diagonal  line  in  Fig.  2  is  due  to  the 
decrease  in  oxygen  concentration  as  a  function  of  altitude.  Because  the 
three-body  attachment  is  proportional  to  the  square  of  the  oxygen 
concentration,  plasma  lifetime  increases  by  two  orders  of  magnitude  for  each 
factor  of  10  decrease  in  atmospheric  pressure. 

As  the  concentration  of  oxygen  diminishes  with  increasing  altitude,  two- 
body  attachment  and  recombination  processes  begin  to  dominate  three-body 
attachment  to  oxygen.  These  processes  include  dissociative  attachment  with 
oxygen  and  ozone,  and  dissociative  recombination  with  positive  ions.  Because 
attachment  depends  on  the  concentration  of  oxygen  and  ozone  but 
recombination  does  not,  recombination  with  positive  ions  is  the  dominant  loss 
mechanism  at  high  altitude.  Consequently,  the  plasma  lifetime  eventually 
becomes  height  independent,  which  results  in  the  horizontal  curves  in  Fig  1. 

For  ng  >  10  cm”  ,  dissociative  attachment  of  electrons  with  positive 
ions  and  then  three-body  electron-electron  attachment  to  positive  ions 
dominate  the  deionization  process.  The  nearly  horizontal  curves  in  Fig.  2  and 


the  progressively  larger  separation  between  curves  indicate  that  three-body 
attachment  (proportional  to  n^)  i3  starting  to  dominate  for  n  >  1014  cm-". 

V  V 

The  curves  in  Fig.  2  were  truncated  at  an  altitude  Tor  which  the  energy 

required  for  ionization  would  increase  the  ambient  air  temperature  by  50*  K. 

12  —3 

For  n  <  10  cm  '  the  amount,  of  beating  below  200,000  ft  is  negligible. 

V 

Noble  Gases:  Significantly  longer  lifetimes  are  possible  with  helium. 

Helium  does  not  form  negative  ions,  so  a  major  electron  attachment 

mechanism  is  minimized.  Because  high  purity  helium  is  expensive  and 

impurities  always  leak  into  a  large  vessel,  helium  with  trace  amounts  of  air  is 

the  most  likely  gas  mixture.  Nitrogen  and  oxygen  impurities  rapidly  form 

positive  ions  in  charge  transfer  collisions  with  helium  ions  and  in  turn 

recimbine  with  electrons  and  so  deionize  the  helium  plasma. 

Simulations  of  plasma  lifetime  in  helium  with  1-ppm,  lO—ppm,  and 

IQO-ppm  air  demonstrate  lifetimes  orders  of  magnitude  longer  than  air  plasma 

1  1 

lifetimes,  as  shown  in  Figs.  3,  4,  and  5,  respectively.  For  ne  2  10  cm 
and  impurity  concentrations  below  100  ppm,  the  lifetime  in  helium  is  ~10 
times  longer  than  an  air  plasma  at  sea  level.  At  an  altitude  of  110,000  ft 
(30.5  km,  8.5  torr),  air  and  helium  plasmas  have  similar  lifetimes.  For  a 
high-purity  helium  plasma,  Deioche  (17)  observed  plasma  lifetimes  longer  than 
those  estimated  in  Fig.  3. 

An  engineering  estimate  for  r  in  a  noble  gas  is  obtained  by  examining 
the  recombination  kinetics  to  determine  if  there  is  a  limiting  reaction.  For 
ne  <  10  cm  ’  in  a  noble  gas  with  a  trace  of  air,  two-body  dissociative 
recombination  is  likely  to  characterize  the  overall  recombination  rate,  and  the 
plasma  lifetime  reduces  to  r  =  l/(kne),  where  k  is  a  two-body  rate  constant. 
Because  dimer  formation  and  charge  transfer  are  fast  compared  to  dissociative 
recombination,  the  dominant  ion  for  recombination  depends  on  the  purity  of 
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tne  noble  gas.  Charge  transfer  resulting  in  N^  and  0^  will  increase  the  rate 
of  recombination  by  a  factor  of  10  compared  to  Heg-  Two-body  rate 
constants  for  noble  gases  and  air  [18]  are  shown  in  Table  III. 


TABLE  III 

DISSOCIATIVE  RECOMBINATION 


Reaction 

3 

Rate  constant  (cm  /s) 

He^  +  e  ->  He  +  He 

O 

X 

O 

pH 

Ne"^  +  e  -*  Ne  +  Ne 

1.8  x  10  7 

Ar^  -f  e  Ar  +  Ar 

9.1  x  1(T7 

Kr+  +  e  -  Kr  +  Kr 

o 

rH 

X 

<£> 

t—i 

Xe^  +  e  “♦  Xe  +  Xe 

2.7  x  Kf6 

N"£  +  e  -*  N  +  N 

3.5  x  10-7 

0+  +  e  -*  0  +  0 

1.9  x  10~7 

Data:  Temperature  =  300*  K,  1  cm3/s  =  10-6  m3/s 


IV.  APPLICATIONS 

Figure  1  sketches  the  use  of  plasma  as  a  reflector  (artificial  ionospheric 
mirror)  and  as  an  absorber.  If  the  high-power  pulses  of  [4]  and  [5]  result  in 
the  formation  of  a  plasma  patch,  will  the  patch  reflect  or  transmit  RF?  If 
RF  breakdown  results  in  a  plasma  with  a  rapid  transition  from  free  space  to 
plasma,  a/X  <  1,  the  plasma  will  reflect.  Or,  if  the  transition  is  diffuse,  the 
piasma  will  reflect  at  a  grazing  angle.  If  the  plasma  gradient  is  too  long,  that 
is,  a  >  A/(4t),  then  it  will  transmit  more  RF  than  it  reflects. 
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As  an  example  of  an  ionospheric  mirror,  consider  a  plasma  at  an  altitude 

of  230,000  ft  (70.1  km),  a  pressure  of  40  mtorr,  an  incident  angle  of  75 

£  _ 1 

degrees  off  broadside,  and  y  =  7.4  *  10  s  .  The  curves  in  Fig.  6  quantify 
the  dependence  of  the  reflection  coefficient  on  a,  the  plasma  to  free-space 
transition,  and  on  frequency.  The  value  of  nQ  =  10  cm  was  chosen  so  as 
to  achieve  a  high  reflection  coefficient  at  a  sharp  boundary  for  frequencies  up 
to  100  MHz.  Note  the  transition  from  a  good  reflector  to  a  poor  reflector  as 
the  gradient  varies  from  sharp,  a  =  10  cm,  to  diffuse,  a  -  1G0  m.  Increasing 
n0  increases  the  maximum  frequency  for  high  reflectivity.  Lowering  the 
altitude  increases  absorption  due  to  collisions,  decreases  plasma  lifetime,  and  so 
requires  more  power. 

As  an  example  of  a  broadband  absorber,  consider  a  plasma  in  air  at  sea 
level  with  v  =  417  *  108  s”1  and  n0<r  =  6  *  1010  cm-3m.  These  values 
were  selected  to  yield  a  maximum  round-trip  absorption  of  40  dB  in  Fig.  7. 
Round-trip  absorption  refers  to  the  sum  of  reflection  from  the  gradient  (3) 
plus  absorption  (8)  due  to  propagating  (h  =  5.888a)  into  the  plasma  and 
reflecting.  Doubling  n  will  double  absorption  (7). 

V 

Absorption  bandwidth  depends  on  the  collision  rate  and  the  plasma 
gradient.  The  high-frequency  cutoff  is  set  by  v  and  has  a  value  of 
The  variation  of  absorption  with  a  in  Fig.  7  indicates  that  high  absorption  at 
a  low  frequency  requires  a  free-space  transition  of  4<r  2  A,  that  is,  one  or  two 
wavelengths.  The  low-frequencv  cutoff  is  ~c/(4a). 

A  plasma  with  a  >  A/(4x)  and  sufficient  electron  number  density  has 
the  characteristics  of  a  broadband  absorber.  This  plasma  absorber  is 
lightweight,  exhibits  high  attenuation  per  wavelength,  has  a  broad  bandwidth, 
and  can  be  switched  on  and  off.  These  characteristics  provide  a  means  to 
control  electromagnetic  scattering.  A  plasma  at  atmospheric  pressure  could  be 
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used  as  an  absorber  for  antenna  pattern  and  b&ckscatter  measurements;  modify 
the  radiation  pattern  of  an  antenna;  act  as  a  high-power  waveguide  switch  or 
crowbar;  and  reduce  edge  diffraction.  Although  these  applications  are 
predicted  and  experimental  evidence  exists  to  verify  the  reflective/absorptive 
characteristics  of  a  collisional  plasma,  an  efficient  means  of  generating  a 
collisional  plasma  is  still  evolving. 


IV.  PRODUCTION  AND  USE  OF  PLASMA 

A.  Plasma  Generation 

Two  requirements  of  plasma  generation  for  an  absorber  are  that  (1)  ng  is 
sufficiently  high  near  the  source,  and  (2)  na  diminishes  as  a  function  of 
distance  away  from  the  source.  A  variety  of  sources  fulfill  these  requirements. 
Electron-beam  impact  ionization  and  photoionization  of  a  seed  molecule  are 
two  of  many  laboratory  techniques. 

A  high-energy  (ICO  keV  to  600  keV)  pulsed  electron  beam  can  generate 
a  plasma  in  air  or  noble  gases.  As  electrons  emerge  from  a  source,  they 
diverge  and  interact  with  matter.  The  spatial  distribution  is  complex  [19]  and 
involves  electron  backscatter,  secondary  electron  generation,  and  multiple 
shallow-angle  scattering.  These  effects  provide  a  mechanism  for  electrons  to 
curve  around  and  fill  in  regions  not  directly  illuminated  by  the  beam. 

Electron  density  from  a  point  source  falls  off  as  1/r  ,  and  electrons  have 
a  maximum  range  due  to  their  interactions  with  matter.  Because  of  this 
continual  interaction,  electron  energy  decreases  with  range.  Below  20  keV,  the 
rate  of  energy  deposition  increases  significantly.  This  increase  results  in  a 


15 


small  increase  in  ng  near  its  maximum  range.  High-energy  electron-beam 
sources  are  commercially  available  and  have  been  used  in  laboratory 
experiments. 

Another  ionization  technique  is  photoionization  of  a  seed  gas  such  as  an 

alkali  or  organic  vapor.  For  efficient  photoionization,  the  medium  must  be 

optically  thick.  To  achieve  the  required  alkali  vapor  density  for  a  1-  to 

4-meter  optical  thickness,  it  must  be  heated  and  contained  in  a  membrane  not 

affected  by  the  heat  or  vapor.  An  alternative  is  to  use  a  readily  ionized 

organic  vr.por  With  a  high  vapor  pressure  such  as  TMAE.  Rewick  [10]  has 

identified  a  number  of  common  materials  that  could  contain  TMAE  vapor. 

There  are  several  commercial  UV  flash  lamps,  rich  below  200  run,  that  could 

readily  photoionize  TMAE  vapor.  A  plasma  generated  via  photoionization 
o 

would  have  a  1/r  distribution  that  is  also  attenuated  because  of  the  optical 
range  of  the  medium. 

B.  Power  Estimation 

An  approximation  for  the  power  per  unit  volume  required  to  sustain  a 
plasma  is  P/V  —  (n  E.)/r,  where  P  is  the  continuous-wave  (CW)  power 

v  1 

(deposited  in  the  plasma),  V  is  the  plasma  volume,  n  is  the  electron-plasma 
density,  E.  is  the  energy  to  generate  an  electron-ion  pair,  and  r  is  the  plasma 
lifetime  (Section  III  B).  For  a  plasma  generated  by  electron-beam  impact,  the 
value  of  Ej  if*  41.5  eV  for  He,  38.2  eV  for  Ne.  26.2  eV  for  Ar,  24.3  eV  for 
Kr,  21.9  eV  for  Xe,  34.6  eV  for  N2,  31.8  eV  for  02,  and  33.7  eV  for  dry  air. 
ll  photoionization  of  a  seed  molecule  like  TMAE  is  the  means  of  ionization, 
then  Ej  is  approximately  7  eV. 
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Consider  two  examples  in  air:  a  plasma,  absorber  at  sea  level  and  an 

ionospheric  mirror.  For  a  plasma  with  n  =  10  cm  '  at  760  torr,  the 

3 

plasma  lifetime  is  1.0  ns,  and  P/V  is  54  MW/m  for  impact  ionization.  An 

ionospheric  mirror  at  230,000  ft  (70.1  km)  with  n  =  10  cm  and 

3 

r  -  800  ms  has  a  power  requirement  of  68  pW/m  delivered  to  the  plasma. 
Plasma  lifetimes  in  a  noble  gas  at  760  torr  an:  orders  of  magnitude  longer 
than  in  air.  which  significantly  reduces  the  power  requirement. 

The  power  per  unit  volume  for  a  noble-gas  plasma  dominated  by 
two-body  recombination  reduces  to 

-  =  k  He  Ej  (9) 

y  e  1 

The  rate  constants  for  noble  gases  are  shown  in  Table  III.  Tbe  rate  constant 
for  dissociative  recombination  of  organic  vapor  ions  like  TMAE+  have  not 
been  measured  but  are  estimated  to  be  of  the  order  of  Xe£  in  Table  III. 

The  power  estimates  in  fable  IV  are  for  impart  ionization  of  helium 
(Ej  2  42  eV)  and  photoionization  of  organic  vapors  (Ej  2  7  eV). 


C.  Use  as  Absoroer 

Even  though  none  of  the  ionisation  sources  mentioned  aboye  generates  an 
Epstein  profile,  tne  a  parameter  in  (2)  can  be  adjusted  tc  approximate  other 
distributions.  The  value  of  a  is  set  to  match  the  5%  and  95%  values  of 
electron  density  produced  by  the  ionization  source.  By  making  this 
adjustment,  an  Epstein  profile  can  be  used  to  estimate  electromagnetic  effects. 
The  alternative  is  to  use  the  exact  electron  distribution  and  a  computational 
model  for  scattering  from  a  plasma  gradient  [?].  An  Epstein  profile  provides  a 
simple  first  estimate. 
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TABLE  IV 

POWER  REQUIRED  TO  SUSTAIN  A  PLASMA 


Plasma 

Density 

-) 


(cm  ^ 


3 

Power  per  Unit  Volume  (W/m  ) 

k  =  10"8  cm3/s  10"7  cm3/s  UT6  cm3/s 

Ej  =  42  eV  7  tV  7  eV 


10B 

70 

X 

10"9 

100 

X 

\Q~° 

1 

X 

10"6 

107 

7 

X 

iO"6 

10 

X 

10"6 

100 

X 

10-6 

108 

700 

X 

10"6 

1 

X 

1Q“3 

10 

X 

IQ-3 

IQ9 

70 

X 

10-3 

100 

X 

10“3 

1 

1C10 

7 

10 

100 

HI11 

700 

1 

X 

103 

10 

X 

103 

1012 

70 

X 

103 

100 

X 

103 

1 

X 

106 

The  data  in  Table  IV  has  been  recast  in  Table  V  to  provide  an  estimate 
of  power  per  unit  area  as  a  function  of  absorption,  bandwidth,  and  gas 
composition.  Bandwidth  was  deduced  from  several  plots  like  Fig.  8,  tor  which 
the  collision  rate  was  reduced  from  that  of  air  at  sea  level  to  that  of  helium 
and  a  eon /argon  at  760  torr.  Parameters  that  yield  predictions  of  attenuation 
up  to  100  dB  are  noted  in  Fig.  8.  For  values  of  absorption  beyond  100  dB 
backscatter  from  local  inhomogeneities  in  the  plasma  gradient  adds  to 
round-trip  absorption  at  low  frequencies.  Predictions  of  round-trip  absorption 
a,  excess  of  100  dB  must  take  these  local  iahomogeneities  into  account. 
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TABLE  V 


POWER  REQUIRED  FOR  VHF  TO  S-BAND  ABSORPTION 


Absorption 

(dB) 

Power 

’  (W/m2) 

Bandwidth  (MHz) 

Helium 

Neon/ Argon 

Helium 

Neon/ Argon 

40 

13,500 

180 

■MRPM 

30 

30,400 

420 

100  -  9,000 

SO 

54,000 

780 

120  -  6,000 

100 

77,780 

1,058 

150  -  6,000 

t  Power  refers  to  CW  power  per  square  meter,  for  a  plasma  with  a  4-m 
thick  plasma  that  approximates  an  Epstein  profile. 


Note  the  trade-off  in  Table  V  between  power  and  bandwidth  for  a  neon/ 

argon-based  plasma  compared  to  helium.  The  high-frequency  cutoff  is 

2 

proportional  to  ir,  power  is  proportional  to  n  ;  and  absorption  (7)  is 
proportional  to  (ne/r).  Because  the  value  of  u  in  neon/argon  is  approximately 
a  factor  of  10  less  than  helium,  the  value  of  n  in  neon/argon  is  a  factor  of 

V 

10  less  than  in  helium  for  the  same  absorption,.  This  reduces  bandwidth  by  a 
factor  of  10,  but  power  decreases  by  a  factor  of  100. 

The  power  in  Tables  IV  ;jcd  V  is  the  power  deposited  in  the  plasma  and 
so  does  not  account  for  source  efficiency.  Generating  plasma  intermittently 
with  a  low  uufc.y  ratio  is  one  way  to  produce  a  large  volume  of  plasma 
absorber  with  simple  equipment.  Other  more  efficient  means  of  ionization 
could  increase  the  duty  ratio  or  permit  the  use  of  gases  that  demonstrate  a 
broader  bandwidth. 
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V.  CONCLUSION 


The  theory,  experimental  evidence,  and  calculations  support  use  of 
plasmas  at  atmospheric  pressure  as  reflective/absorptive  media. 

Electromagnetic  characteristics  were  evaluated  with  an  Epstein  profile  model 
that  can  be  used  to  approximate  the  electron  distributions  generated  by  an 
ionization  source.  Although  the  modeling  of  the  source  distribution  function  is 
approximate,  the  Epstein  profile  model  provides  a  simple  estimate  of 
Electromagnetic  absorption  and  scatterir  g. 

Plasma  generated  in  air  or  heliun  at  atmospheric  pressure  has  the  novel 
property  of  being  an  excellent  broadb«jid  absorber  from  VHF  to  X-band.  The 
power  required  to  sustain  a  plasma  ii  a  loble  gas  like  helium  is  orders  of 
magnitude  less  than  >n  air,  because  aob) ;  gases  do  not  form  negative  ions.  A 
neoii/argon  plasma  has  the  property  of  nigh  absorption  from  VHF  to  S-band 
an  piircs  less  power  to  sust  in  that  a  helium  plasma,  because  its 
momentum-transfer  collision  r  ■>  ce  i  i  lover. 

^hese  characteristics  suggest  that  plasmas  at  atmospheric  pressure  have 
the  potential  of  becoming  a  useful  engineering  material  from  VHF  to  X-band. 
Ine  gap  in  knowledge  to  thii  practical  end  involves  a  better  understanding  of 
scattering  parameters,  ionL:-  .ion  and  recombination  gas  kinetic:,  and  efficient 
ionization  techniques. 


ACKNOWLEDGMENT 

This  work  was  supported  by  Air  Force  Office  of  Scientific  Research, 
AFQSR/NP,  contract  F49620-85-K-0013.  The  author  acknowledges  many 
productive  discussions  with  D.  Eckstrom  and  W.  Che  nut. 


20 


REFERENCES 


[1]  N.  C.  Gerson,  Radio  Wave  Absorption  in  the  ionosphere.  New  York: 
Pergarnon  Press,  1962. 

[2]  M.  Gunar  and  R.  Mennella,  "Signature  Studies  for  a  Re-Entry  System," 
Proceeding?  of  the  Second  Space  Congress— New  Dimensions  in  Space 
Technology.  Canaveral  Council  of  Technical  Societies,  pp.  515-548,  1965. 

[3]  S.  Glasstone  and  P.  J.  Dolan,  The  Effects  of  Nuclear  Weapons. 
Washington:  United  States  Department  of  Defense  and  Energy  Research 
and  Development  Administration,  pp.  461-513,  1977. 

[4]  N.  D.  Borisov  and  A.  V.  Gurevich,  "High-Frequency  Pulsed  Air 
Breakdown  in  Intersecting  Radio  Beams,"  Geomagnetism  and  Aeronomv. 
vol.  20,  pp.  587-591,  1980. 

[5]  A.  V.  Gurevich,  "An  Ionized  Layer  in  a  Gas  (in  the  Atmosphere),"  Sov. 
Phvs.  Usp..  vol.  23,  pp.  862-865,  1980. 

[6]  P.  S.  Epstein,  "Reflection  of  V.  ^ves  in  an  Inhomogeneous  Absorbing 
Medium,"  Proc.  Nat.  4cad  Sci..  Wash,  vol.  16,  pp.  627-637,  1930. 

[7j  K.  G.  Budden,  The  Propagation  of  Radio  Waves.  Cambridge  (Great 
Britain).  Cambridge  University  Press,  pp.  470-475,  550-582,  1985. 


21 


I 


[8]  D.  F.  Anderson,  "A  Photoionization  Detector  for  the  Detection  of  Xenon 
Light/'  IEEE  Transactions  on  Nuclear  Sciences,  voi.  NS-28, 

pp.  842-848,  1981. 

[9]  R.  A.  Holroyd,  J.  M  Treses,  C.  u.  Woody,  and  R.  A.  Johnson, 
"Measurement  of  the  /Absorption  Length  and  Absolute  Quantum 
Efficiency  of  TMAE  aud  TEA  from  Threshold  to  120  nm,"  Nucl. 

Instrum.  Methods  PKs.  Res.  A  ( Netherlands!,  vol.  A261,  pp.  440-444, 

1987. 

flOj  R.  T.  Rewick,  M.  L.  Schumacher,  S.  L.  Shapiro,  T.  B.  Weber,  and  M. 
Cavalli-Sforza,  "Tetrakis{dimethylamino)etiiyiene:  Identification  of 
Impurities  and  Compatibility  with  Common  Metal,  Polymer,  and  Ceramic 
Laboratory  Materials,"  Analytical  Chemistry,  vol.  60,  pp.  2095-2099, 

1988. 

(Mj  M.  Mitchner  and  C.  Kruger,  Partially  Ionized  Gases.  New  York:  Wiley 
and  Sons,  pp.  47-53,  155-162,  1973. 

[12]  B.  Tanenbaum,  Plasma  Physics.  New  York:  McGraw-Hill,  pp.  62-85, 

1967. 

[13]  C.  Lanczos,  "A  Precision  Approximation  of  the  Gamma  Function," 

L  SIAM  Numerical  Analysis.  Ser.  B,  vol.  1,  pp.  85-96,  1964. 


22 


[14]  Y.  Itikawa,  "Effective  Collision  Frequency  of  Electrons  in  Gases," 

(Note,  typo  in  Table  II,  10“18  should  be  Iff8  secrlcm\  The  Physics  of 
Fluids,  voi.  16,  pp.  831-835,  1973. 

[15]  R.  J.  Vidmar,  "Generation  of  Tenuous  Plasma  Clouds  in  the  Earth's 
Atmosphere,"  Annual  Report  for  AFOSR/NP  contract  F49620-85-K-0G 13, 
SRI  International,  Menlo  Park,  CA,  pp.  43,  1987. 

[16]  R.  J.  Vidmar,  "Generation  and  Properties  of  Tenuous  Plasma  Bodies  at 
Atmospheric  Pressure,"  SRI  Annual  Report  for  AFOSR/NP  contract 
F49620-85-K-0013,  SRI  International,  Menlo  Park,  CA,  pp.  35,  1988. 

[17]  R.  Delcche,  R.  P.  Monchicourt,  M.  Cheret,  and  F.  Lambert, 
"High-Pressure  Helium  Afterglow  at  Room  Temperature,"  Physical 
Review  A.  vol.  13,  pp.  1140-1176,  1976. 

[18]  l.  G.  Christophorou,  ElsarpR-Moleculs  InterartigM_and  Their 
Applications.  Volume  2.  Orlando:  Academic  Press,  pp.  74-77,  1984. 

[19]  F.  R.  Attix,  Introduction  to  Radiological  Physics  and  Radiation 
Dosimetry.  New  York:  Wiiey,  pp.  160-187,  1986. 


23 


FIGURE  CAPTIONS 


Fig.  1.  Plasma  at  atmospheric  pressure  used  (a)  as  a  reflector  and 
(b)  as  an  absorber. 

Fig.  2.  Plasma  lifetime  of  air  versus  altitude. 

Fig.  3.  Plasma  lifetime  of  helium  with  1  ppm  air. 

Fig.  4.  Plasma  lifetime  of  helium  with  10  ppm  air. 

Fig.  5.  Plasma  lifetime  of  helium  with  100  ppm  air. 

Fig.  6.  Reflection  from  an  Epstein  profile  at  230,000  ft  (70.1  km). 

Fig.  7.  Absorption  plus  reflection  from  an  Epstein  profile  in  air  at  760  torr. 

Fig.  8.  Absorption  plus  reflection  from  an  Epstein  profile  in  helium  and 

neon/argon  at  760  torr. 
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FIGURE  1  PLASMA  AT  ATMOSPHERIC  PRESSURE  USED 
(8)  AS  A  REFLECTOR  AND  (b)  AS  AN  ABSORBER 
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FIGURE  2  PLASMA  LIFETIME  OF  AIR  VERSUS  ALTITUDE 
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FIGURE  5  PLASMA  LIFETIME  OF  HELIUM  WITH 
100  ppm  AIR 
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FIGURE  6  REFLECTION  FROM  \N  tPSTEIK  pRO”ILE  AT 
230  kft  (70  km) 
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(b)  NEON/ ARGON 


FIGURE  8  ABSORPTION  PLUS  REFLECTION  FROM  AN 
EPSTEIN  PROFILE  IN  HELIUM  AND 
NEON/ARGON  AT  760  TORR 
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ABSTRACT 

A  streseed-skin  inflatable  target  support 
provides  an  improvement  over  a  foam  column 
for  radar  crosa  section  (RCS)  seaaureaenta 
in  an  anechoic  chamber .  Theoretical 
analysis  indicates  that  backscatter  from 
the  support  is  minimi?.  »d  because  its  aass 
is  reduced  beJov  that  of  a  foaa  coluan-and 
is  distributed  co  favor  incoherent 
scattering.  Compared  with  a  foaa  column, 
a  pressurized  thin  shell  has  superior 
mechanical  stability  under  both  axial  and 
transverse  loads.  Experimental 
observations  using  Mylar — a  low  dieisctric 
constant,  high  tensile  strength  fila — 
confirm  these  results.  Spurious 
reflsctions  froa  rotational  aachinsry 
located  below  an  inflatable  column  are 
reduced  by  a  layer  of  absorber  within  the 
base  of  the  inflatable  support. 

Keywords:  RCS  Measurement,  Target  Support, 
Stressed  Skin,  Inflatable  Column,  Cluttsr 
Reduction 


1 .  INTRODUCTION 

Anechoic  chamber  clutter  1 iaits  both 
precision  and  noiss  floor,  making  it 
difficult  for  RCS  engineers  to  sake 
accurate,  low-revel  RCS  measurements.  A 
principal  sourca  of  cluttsr  is  the  target 
support  column.  A  target  support  column 
must  meet  both  electromagnetic  and 
mechanical  requirements.  It  must  have  low 
RCS  over  a  br<.  ad  bandwidth  and  stably 
support  a  prescribed  weight.  For  accurate 
measurement  of  vector  quantities,  it  is 
important  that  the  target  support  be  able 
to  return  to  a  reference  position  a'ter 
movement . 


Traditional  target  support  methods  include 
foam  columns,  ovei  sad  strings,  and  mstal 
pylons  [1].  A  fo<  .  column  has  difficulty 
supporting  a  heavy  load  while  maintaining 
a  low  backscatter,  and  viscoelastic 
properties  affect  its  positional  accuracy. 
A  string  support  system  is  difficult  to 
control  and  requires  an  attachment  point 
to  the  target.  Met *1  ogival  pedestals  can 
support  heavy  loads,  but  have  a  larger  RCS 
than  other  support  methods. 

An  inflatable  column  offers  both  low  RCS 
and  heavy  lift  capacity.  An  inflatable 
column  was  designed  at  RATS CAT  in  1964  [2} 
using  nylon-coated  neoprene.  The  column, 
designed  to  support  200  lb  in  a  40-knot 
wind,  was  inflatsd  to  0.2  psi  and  was 
stable  enough  to  permit  vector 
measurements .  Performance  of  the  RATS CAT 
column  was  comparable  to  that  of  foaa. 
Current,  advances  in  aatsrial  science  have 
given  us  the  ability  to  fabricate 
inflatable  structures  with  properties  that 
•re  superior  to  foaa  construction. 

An  inflatable  column  has  several 
advantages  over  a  foam  column  of  identical 
■hape.  First,  an  inflatabla  column  has  a 
lowsr  RCS  than  a  foaa  column.  Second,  an 
inflatable  column  made  of  a  membrane 
material  that  is  significantly  loss 
viscoelastic  than  foam  can  havs  superior 
mechanical  properties,  which  results  in 
the  ability  to  conduct  angular 
measurements  in  a  sore  stable  manner  and 
to  support  heavier  loads.  Finally, 
ground-related  clutter  can  be  reduced  over 
a  broad  bandwidth  by  filling  the  base  of 
the  column  with  broadband  absorber. 

Section  2  addresses  design  parameters  for 
an  inflatable  column,  including  both 
electrical  and  mechanical  considerations. 
Section  3  describes  fabrication  details 
that  resulted  in  a  prototype  structure. 
Section  4  discusses  anechoic  chamber 
Measurements ,  ) comparing  the  scattering 
froa  a  square  plate  mounted  at  45*  on  a 
foam  column  with  scattering  from  the  aame 
plate  mounted  atop  the  inflatable  column. 
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U.-5SIGK  CONSID2RATTGHS 

Ti.s  tio  principal  requirements  fc  r  a 
support  ce.luan  era  Minimum  RC5  ova*  a 
specified  frequency  rang®  and  stotble 
mechanical  support:  up  to  a  specified 
weight..  Tesign  parawat&rs  fa::  fofe*  and 
inflatable  columns  includes 
eluotroBa^nstlc  scattering,  axial  loading, 
kr*c*/er*c  loading,,  and  torsional 
rigidity.  An  inflatable  column  can  be 
deeign.vd  with  minimum  RCS  and  tfc«  ua me  or 
tuperi-v..  ae.il' ntical  performance  compared 
with  a  foam  eo^an. 


Figure  l  is  a  plot  of  scattering  per  unit 
length  as  a  function  of  frequency  for  a 
12-in.  diameter  foam  cylinder  of  • 
permittivity  1.04  (dashed  curve)  and  a  12- 
in.  diameter  Mylar  shell  of  thickness 
0.002  in.  and  permittivity  2.8  (solid 
curve) .  At  low  frequencies,  scattering 
from  the  inflatable  structure  is 
considerably  reduced  coapared  with  that 
froa  a  foaa  column.  At  sufficiently  high 
frequencies,  this  inflatable  column  would 
exhibit  a  high  RCS  because  the  film 
thickness  approaches  a  significant 
fraction  of  a  wavelength. 


2 . 1 .  Elec  roaegis^t J,c  Properties 

The  electromagnet  to  scattering  free  a 
cesluan  can  be  *sti*iat«d  by  considering  the 
cwo-dimcnsiohii.;.  scattering  free  &  class  of 
infinite  cylinders  j[  3  *i .  A  model  for 
scattering  fro*  *  fee®,  column  is  a  solid 
dielectric  tylin&sr  whose  permittivity  is 
nearly  unity.  The  scattering  cross 
section  per  .mit  length  is: 
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where  a  is  the  radius  of  tie  cylinder, 
ko  is  the  Jfree  tpece  rave  number,  and 
er  is  the  relative  pmrrittivity  of  the 
cylinder,  T h*  response  .  z  the  .uum  for 
either  parallel  cr  pwEpcnc iculai 
polarisation. 


For  n  inflatable  coImish,  similar 
expressions  can  toe  derived,  whore  the 
model  is  noj  i  cy" indrical  shell.  The 
scatter :  nq  cruse  section  per  unit,  length 
is: 
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where  c„  -  1  when  n  -  0  and  »,  *>  2  when 
n  >  0,  and  d  Is  the  shell  thickness. 
Slightly  different  values  ar«  obtwir.ed  '’or 
parallel  and  perpendicular  polariretion. 


Fig.  1  Normatred  BacfcscaBer  Crow  Section  par  Unit  len^h  tiww 
iriinrtely  Long  CylncMw* 

cone-shaping  of  either  column  can  be  usod 
to  further  reduce  the  direct  backscatter . 

2.2  Mechanical  Properties 

The  model  for  comparison  of  mschanical 
propertise  is  a  column  that  is  fixed  at 
one  end  and  free  to  move  at  the  other. 

The  maximum  buckling  load,  W,  for  a  foam 
column  is  givsn  by  (4): 

w  - 

“I  4CT  * 


where  a  .!«  the  radius  of  tha  column,  L  is 
the  height,  and  E  is  the  alaetic  modulus. 
Ths  maximum  load  for  an  inflatable  column 
is  determined  by  the  cross-sectional  area 
at  the  top  and  the  pressure  in  the  column 


* 


The  pressure,  P,  in  the  column  is  limited 
by  ths  maximum  hoop  stress,  S  -  Pa/d, 
(which  occurs  at  maximum  radius) .  The 
maximum  load  Is:  X  •  xadS,  where  d  is  the 
film  thicicness.  For  taller  columns,  it  is 
appropriate  to  consider  failure  due  to 
buckling: 

w  -  ^4^ 

4tr 

A  10- ft  foam  column  of  6- in.  radius  and 
elastic  modulus  of  1000  psi  is  capable  of 
supporting  43  io«,  with  a  safety  factor  of 
4.  A  10-ft  air  column  of  6-in.  radius 
would  need  to  be  inflated  to  0.76  pci  to 
support  a  similar  load.  For  a  cylindrical 
column,  this  pressure  would  exert  a  4.5 
lb/ in.  hoop  stress  (Sd)  at  the  base.  A 
safety  factor  of  two  would  require  the 
cylinder  to  withstand  9  Ih/in.  For  a 
conical  column  with  a  bnse  r  idius  of  16 
in. ,  the  stress  at  the  base  would  require 
a  design  of  12  lb/in.  This  can  he 
accomplished  using  2-mil  Mylar.  (For 
comparison,,  about  20  mils  of  polyethylene 
film  would  be  needed  to  me*t  the  same 
stress  requirement.) 

The  deflection  of  a  column  under 
transverse  loading  can  be  determined  to 
first  order  using  email  deflection  theory. 
For  a  conical  foam  column,  the  maximum 
deflection  is: 


Although  the  foam  and  inflatable  columns 
exhibit  similar  characteristics,  foam  is 
much  more  viscoelastic  than  Mylar. 
Consequently,  a  full  comparison  between 
foam  and  inflatable  construction  must 
treat  how  the  elastic  modulus  varies  as  a 
function  of  time,  temperature,  and  stress. 
For  example,  the  elastic  modulus  may 
decrease  by  a  factor  of  10  in  1  hour  while 
a  high  stress  is  applied.  This  variation 
in  E  with  time  results  in  a  time-dependent 
deflection.  Foam  may  exhibit  an 
additional  deflection  of  10%,  while  Mylar 
creep  is  an  order  of  magnitude  leas. 

Under  severe  loading,  a  foam  column  may 
nevar  return  to  its  reference  position 
after  the  transverse  load  is  removed. 

The  response  of  a  cylindrical  column  under 
torsion  can  be  analyzed  using  the 
following  equation  [5]: 

..  TL(!  +  v) 

»  **  — ^BT — 

where  e  is  the  angular  displacement,  T 
is  the  applied  torque,  I  is  the  moment  of 
inertia  for  a  body  of  revolution  (Ic  and 
I.) ,  and  v  is  Poisson's  ratio.  A 
comparison  of  the  two  columns  indicates 
that  the  foam  column  exhibits  an  angular 
displacement  of  0.0016°/ft~lb#  while 
the  inflatable  column  changes 
0 . 002 4®/ ft -lb. 
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where  Ar  »  i2  -  r, ,  r,  is  the  radius 
at  the  top,  r?  is  the  radio  at.  the  base, 
Ic  -  fr  4  '*  for  a  solid  circle,  ar.d  F  is 
the  tr-  .  ”  ie  force. 


Although  the  deflections  and  angular 
displacements  of  this  version  of  the 
inflatable  column  are  slightly  higher  than 
those  of  a  foam  column,  the  inflatable 
returns  to  its  reference  position, 
reliably.  The  reliability  of  return  opens 
up  the  possibility  of  high  accuracy 
seasurements  involving  vector  subtraction. 


For  an  inflatable  column,  the  deflection 

is : 
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where  I,  -  sir'd  for  a  thin  annulus. 

A  comparison  of  a  foam  and  inflatable 
construction  using  the  above  dimensions 
indicate*  that  a  10-1 L  transverva  loud 
applied  to  a  ^7-in.  column  will  produce  a 
0.112-in.  deflection  cf  a  foam  column  and 
a  0.332-in.  deflection  of  an  inflatable 
column.  (The  elastic  modulus  of  Mylar  i 
•5 00 , 00  v  psi  ) 


3 .  CONSTRUCTION 

The  cone-shaped  column  in  Figure  2  was 
constructed  using  2-mil  Mylar  film.  The 
column  is  inflated  to  0.75  psi,  using 
nitrogen  gas,  and  can  support  a  load  of  50 
lb.  The  vertical  bag  seams  are  double  lap 
joints  using  an  adhesive  tap*  on  both  the 
inside  and  outside.  Stress  testing 
indicates  that  the  material  will  fail 
before  the  adhesive  joint.  A  heat- 
shrinkable  film  is  applied  to  a  foam  plug 
and  inserted  at  the  top  of  the  column. 

The  internal  pressure  of  the  bag  forces 
the  tapered  foam  plug  to  seal  the  top  of 
the  bag.  Tha  base  is  constructed  using  an 
aluminum- foam  sandwich,  vith  the  foam 
cantilevered  beyond  the  aluminum  plates. 
This  arrangement  peefts  absorber,  placed 
inside  the  column,  to  cover  the  aluminum 
plate.  A  ring  of  flat  laainate  absorber 
is  placed  around  the  outside  of  the  column 
to  complete  the  Ay  seal  to  the  base.  The 
Mylar  ban  is  attached  to  the  base  using  a 
special  adhesive  technique. 


MATERIALS: 

Alummum 

Foam 

Mytetf 

l.aminale  Atucvtvar 
Pyramid  Ateorocsr 


DESIGN  DATA: 

Deity':  Weight:  SO  ii 
Prewura:  C.7S  psi 
Haigm:  97  in. 

To*?:  13  in.  OD 
Base:  32  In.  OO 
Taper:  6.5* 
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Fig  2.  Stres*fl<j-Skin  inflatable  Structure 
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Fig.  3.  Scattering  From  a  45*  Rate 


Two  prototype  column*  have  been 
constructed  end  Inflated  to  the  design 
pressure  of  0.75  p*i.  They  have  operated 
reliably  for  three  Months  and  are 
continuing  to  operate  at  the  tiae  of  this 
submission.  Design  lifetime  is  one  year 
or  sore. 

4 .  PERFORMANCE 

Preliminary  measurements  on  the  cone- 
shaped  inflatable  column  have  been  made  in 
our  anechoic  chamber  [6]  using  an  HP-8510 
automatic  network  analyz  r.  A  coapariaon 
is  made  between  the  foam  column  and  the 
inflatable  column. 

A  45°  plate,  14  in.  on  a  side  is  mounted 
on  the  column,  in  a  way  to  cast  a  direct 
reflection  toward  the  base.  Figure  3  ie  a 
range  domain  plot  of  the  foam  column  and 
inflatable  column  responses.  The 
horizontal  axis  represents  range  in 
centimeters,  the  vertical  axis  is  rsletivs 
amplitude  in  dB,  and  the  frequency  of 
operation  is  600-4600  MHz .  The 
measurement  was  made  using  the  bare  column 
as  a  reference  and  then  subtracting  this 
reference  from  the  measurement  obtained 
with  the  plate  on  ton  of  the  column. 


The  second  peak  is  a  triangular  path 
reprasenting  scattering  from  the  plate  to 
the  baas  of  the  column  and  back  to  the 
horn  (or  vice  verse.)  The  response  occurs 
earlier  in  tiae  for  the  inflatable  column 
because  its  base  is  higher  than  the  base 
of  the  foam  column.  The  amplitude  of  the 
response  for  the  inflatable  column  is 
reduced  by  20  dB  over  the  foam  column 
measurement.  This  reduction  is  due  to  the 
absorber  in  the  base  of  the  air  column. 

The  third  peak  is  determined  by  scattering 
from  the  plate  to  the  base  of  the  column, 
back  to  the  plats,  and  then  returning  to 
the  receive  horn.  For  the  foam  column, 
this  double  bounce  path  is  the  dominant 
scatterar.  For  the  inflatable  column,  the 
double-bounce  response  occurs  earlier  in 
tiae  and  is  35  dB  below  the  foam  column 
result.  Additional  attenuation  of  this 
double  bounce  can  be  obtained  by  using 
better  (taller)  absorber. 

Looking  further  downrange,  multiple 
reflections  from  the  base  of  the  foam 
column  persist  for  soma  time.  This  is  not 
apparent  when  the  inflatable  column  is 
used,  resulting  in  a  reduction  of  about  10 
dB  in  ground-related  clutter. 


Both  curvea  display  three  mein  scattering 
peeks.  Proceeding  downrange  from  the 
illuminating  home,,  the  first  peak 
represents  the  scattering  free  the  plate. 
The  fine  structure  is  scattering  from  the 
tip,  corners,  and  tail.  Small  differences 
between  the  two  curves  is  attributed  to 
placeeent  of  the  45”  plate. 


Although  an  inflatable  column  offers  e 
means  to  reduce  scatter  from  its  base, 
beckscatter  from  the  tep  plug  dominates 
backscattar  from  the  column.  Because  the 
foaa  density  and  dielectric  constant  of 
the  plug  of  the  inflatable  column  are 
similar  to  thoee  of  a  good  quality  foam 
column,  backscattar  pei  unit  length  fro* 
the  plug  is  the  same  as  that  from  *  foaa 
column. 
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Total  backscatter  from  the  inflatable 
column  may  be  significantly  lass  than  or 
equal  to  backscatter  fro*  a  foa*  column, 
depending  on  how  such  of  the  column  is 
illuminated.  If  air  of  the  inflatable 
column  is  illuminated,  backscatter  from 
the  plug  will  dominate  scatter  froa  the 
Mylar  skin.  However,  the  length  of  the 
plug  is  typically  less  than  lOt  of  the 
total  length  of  a  foaa  column.  Hence, 
backscatter  from  a  completely  illuminated 
inflatable  column  would  be  at  least  10  dS 
less  than  that  froa  a  foam  column. 

At  some  facilities  only  a  portion  of  tne 
support  column  is  illuminated.  For 
example,  the  preliminary  measurements 
reported  in  this  paper  were  obtained  in  an 
anechoic  chamber  with  high-gain  horns 
pointed  toward  a  spot  2  ft  above  the 
support  column.  Test  bodies  are  elevated 
to  this  spot  with  foam  fixtures  that  mount 
to  the  support  column.  Backscattsr  was 
primarily  from  the  test  body,  its  fomr 
mount,  and  foam  from  the  top  of  tha 
support  column.  Because  of  this  mounting 
arrangement,  the  amount  of  foam 
illuminated  is  nearly  the  same  for  the 
inflatable  and  foam  columns.  Measured 
backscatter  from  the  two  columns  was 
roughly  equal.  An  improved  version  of  the 
inflatable  column  described  above  would 
include  a  plug  contrived  to  minimize 
backscatter  without  compromising 
structural  integrity. 


column,  the  transverse®  and  torsional 
rigidity  of  th*  column  is  baaed  on  the 
mechanical  properties  of  a  pressurised 
thin-* hell  Mylar  con*.  Mylar  is  much  leas 
viscoelastic  than  typical  foams,  ao 
variations  in  elastic  modulus  with  time, 
temperature,  and  applied  load  are  less  for 
Mylar.  Because  transverse  and  torsional 
deflections  are  inversely  proportional  to 
tha  elastic  modulus,  the  eferenen 
position  of  an  inflatable  column  is  better 
defined  and  can  ba  recovered  after 
transverse  or  torsional  deflections. 

These  electr •'  *1  and  mechanic*!  benefits 
suggest  that,  inflatable  supports  are 
beneficial  to  high-accuracy  vector 
measurements . 
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